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ABSTRACT 



We present an improved analysis of halo substructure traced by RR Lyrae 
stars in the SDSS stripe 82 region. With the addition of SDSS-II data, a revised 
selection method based on new ugriz light curve templates results in a sample of 
483 RR Lyrae stars that is essentially free of contamination. The main result from 
our first study persists: the spatial distribution of halo stars at galactocentric 
distances 5-100 kpc is highly inhomogeneous. At least 20% of halo stars within 
30 kpc from the Galactic center can be statistically associated with substructure. 
We present strong direct evidence, based on both RR Lyrae stars and main 
sequence stars, that the halo stellar number density profile significantly steepens 
beyond a Galactocentric distance of ~30 kpc, and a larger fraction of the stars 
are associated with substructure. By using a novel method that simultaneously 
combines data for RR Lyrae and main sequence stars, and using photometric 
metallicity estimates for main sequence stars derived from deep co-added u-band 
data, we measure the metallicity of the Sagittarius dSph tidal stream (trailing 
arm) towards R.A.~ 2'^-3'' and Dec ~ 0° to be 0.3 dex higher {[Fe/H] = -1.2) 
than that of surrounding halo field stars. Together with a similar result for 
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another major halo substructure, the Monoceros stream, these results support 
theoretical predictions that an early forming, smooth inner halo, is metal poor 
compared to high surface brightness material that have been accreted onto a 
later-forming outer halo. The mean metallicity of stars in the outer halo that 
are not associated with detectable clumps may still be more metal-poor than the 
bulk of inner-halo stars, as has been argued from other data sets. 

Subject headings: methods: data analysis — stars: statistics, RR Lyrae — 
Galaxy: halo, stellar content, structure 



1. Introduction 



Studies of the Galactic halo can help constrain the formation history of the Milky 
Way and the galaxy formation process in general. For example, within the framework 
of hierarchical galaxy formation ( Freeman fc Bland-HawthornI I2OO2I ). the spheroidal com- 



(Johnston, Hernauist & Boltelll996: 


Helmi & White 


1999; 


Bullock. Kravtsov & Weinberg 


2001 



Harding et al.ll200ll ). The number of these substructures, due to mergers and accretion over 
the Galaxy's lifeti me, may provide a crucial test for p roposed solutions to the "missing 
satellite" problem (IBuUock. Kravtsov fc Weinberg||200ll ). Substructures are expected to be 
ubiquitous in the outer halo (galactocentric distance > 15 — 20 kp c), where the dynamical 
timescales are sufficiently long for them to remain spatially coherent (IJohnston. Hernauist fc Bolte 
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2002 ) , and indeed many have been discovered (e.g.. Ilbata et al 
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IVivas fc Zinn 



20081 ) . Various luminous trac- 



ers, such as main sequence turn-off stars, RR Lyrae variables, or red giants are used to detect 
halo substructures, and of them, RR Lyrae stars have proven to be especially useful. 



RR Lyrae stars represent a fair sample of the old halo populati on (ISmithlll995h. They 
arenearly standard candles {{My) = 0.59 ± 0.03 at [Fe/H] = — 1.5. ICacciari fc Clementini 
2OO3I ). are sufficiently bright to be detected at large distances (5 — 120 kpc for 14 < V < 21), 
and are sufficiently numerous to trace the halo substructure with good spatial resolution 
(~ 1.5 kpc at 30 kpc). Fairly compl ete and relatively clean samples of RR Lyrae stars can 
be selected using single-epoch colors ( Ivezic et al. 2005[). and if multi-epoch d a ta existj using 



variabiht v (llvezic et all I2OOOI : QUEST 



SEKBO, Keller et al 



20081 : LONEOS-I, 



Vivas et al.l 



Miceli et al. 



200 ll: ISesar et al.l 120071 . lOe Led I2OO8 
2008h . A useful comparison of recent 



R R Lyrae surveys in t erms of their sky coverage, distance limits, and sample size is presented 
by iKeller et al.l (120081 ) (see their Table 1). Compared to other surveys, the so-called stripe 
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82 survey (an equatorial strip defined by declination limits of ±1.27° and extending fr om 
R.A.~20^ to R.A.~4^) based on Sloan Digital Sky Survey data (SDSS; IVork et alj[200o[ ) is 
the deepest one yet obtained (probing distances to ~ 110 kpc), and the only survey with 
available 5-band photometry. 

In an initial study based on SDSS stripe 82 data, we discovered a com plex spatial 
distri bution of halo RR Lyrae stars at distances ranging from 5 kpc to ~100 kpc (jSesar et al. 



20071 . hereafter S07). The candidate RR Lyrae sample was selected using colors and low- 
order light curve statistics. In this paper, we present a full light curve analysis enabled by 
the addition of new stripe 82 observations. The resulting sample of RR Lyrae stars is highly 
complete and essentially free of contamination. The extended observations now allow the 
determination of flux-averaged magnitudes, thus providing better luminosity and distance 
estimates. This reduces the uncertainty in the spatial distribution of RR Lyrae stars, and 
increases the overall accuracy of computed density maps. 

While this paper was in preparation, another stu dy of RR Lyrae sta rs in stripe 82, 
based on essentially the same data set, was announced (IWatkins et al.ll2009l ). In addition to 
analysis that overlaps with their work, here we also discuss 



• An RR Lyrae template construction for the SDSS bandpasses, and analysis of the 
template behavior 

• Template fitting and visual inspection of best-fits for single-band and color light curves, 
which results in an exceedingly clean final sample 

• Extended color range for the initial selection, which results in about 20% more stars 
than in the sample analyzed by Watkins et al. 

• A step-by-step analysis of the sample incompleteness at the faint end 

• A new method for constraining the metallicity of spatially coherent structures that are 
detected using main sequence and RR Lyrae stars 

• A comparison of the observed RR Lyrae spatial distribution with predictions based on 
an oblate halo model constrained by SDSS observations of main sequence stars. 



We emphasize, however, that both studies obtain consistent main results: strong substruc- 
ture in the spatial distribution of halo RR Lyrae stars and evidence for a steepening of the 
best-fit power-law description beyond a galactocentric radius of ~30 kpc. 

Our data set and the initial selection of candidate RR Lyrae stars are described in § [21 
while the construction of light curve templates is presented in § [31 Various properties of 
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the resulting RR Lyrae sample, with emphasis on the spatial distribution, are analyzed in 
§111 In § [5] we compare the spatial distribution of RR Lyrae stars and main sequence stars 
and constrain the metallicity distribution of the Sgr tidal stream. Our main results are 
summarized and discussed in § O 



The Data 



2.1. An Overview of the SDSS Imaging and Stripe 82 Data 



The Sloan Digital Sky Survey (SDSS; lAdelman-McCarthy et al.ll2008l ) provides homo- 
geneous and deep (r < 22.5) photometry in five bands {u, g, r, i, and z) of more than 7,500 
deg^ of the North Galactic Cap, and three stripes in the South Galactic Cap totaling 740 
deg2. The central stripe in the South Galactic Cap, stripe 82 (20'^32'" < ajsooo.o < 04'^00'", 
— 1.26° < (5j2ooo.o < +1-26°, ~ 280 deg^) was observed repeatedly (not all imaging runs cover 
the entire region) to enable a deep co-addition of the data and to enable the discovery of 
variable objects. It was the largest source of multi-epoch data during the first phase of the 
SDSS (SDSS-I), with a median of 10-20 observations per source (S07, iBramich et al.ll2008l ) 
and it remai ns so thanks to addit ional imaging runs obtained during the SDSS-II Supernova 
(SN) survey ( Frieman et al. I [20081 ). Some of the SDSS-II SN imaging runs were taken in bad 
photometric conditions (e.g., with fast-changing cloud c overage), rendering the photometric 
data unreliable if not recalibrated. We have followed the llvezic et al.l (120071 ) prescription (see 
their Section 3) and recalibrated 180 SDSS-II imaging runs to 1% photometry (magnitude 
zero-point uncertainty < 0.01 mag). The repeated photometric measurements demonstrate 
that the photometric precision is about 0.02 mag at the bright end and reaches 0.05 mag 
at magnitudes of 20, 21.5, 21.5, 21.5, and 19 for the ugriz bands, respectively (S07). The 
morphological inform ation from the images allows reliable star- galaxy separation to r ~ 21.5 
Jhupton et al.ll2002h . 



We have combined SDSS-I imaging runs used in S07 with recalibrated SDSS-II SN 
imaging runs, and created a catalog with ~ 4.4 million unique unresolved (point source) 
objects. Each object in the catalog was observed at least 4 times in the ugr bands, and the 
cadence of re-observations is shown in Figure [H On average, the SDSS stripe 82 objects 
were most often re-observed every two days (the SDSS-II SN Survey cadence), followed 
by 5-day, 10-day and yearly re-observations. We time-average repeated ugriz observations 
(median PSF magnitudes) and corr ect them for interstellar extinction using the map from 
Schlegel. Finkbeiner &: David (119981 ). We also compute various low-order statistics, such 



as the root-mean-square (rms) scatter (a), P^r degree of freedom (x^), and light-curve 
skewness (7) for each ugriz band and each object (see S07 for details). For objects brighter 
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than g = 21, the median number of observations per band is 30. 

The constructed catalog is a superb dataset for the selection and identification of RR 
Lyrae stars. The three-fold increase in the median number of observations per object over 
the S07 catalog allows a reliable measurement of the period of pulsation and enables the 
construction of period-folded ugriz light curves. As we describe in the next Section, RR 
Lyrae stars can be unambiguously identified using their period-folded ugriz light curves, 
and an essentially complete sample of stripe 82 RR Lyrae stars can then be constructed. 



2.2. The Selection of RR Lyrae Stars 

In this section we describe the selection of 483 RR Lyrae stars in SDSS stripe 82. The 
positions and photometry (MJD, ugriz magnitudes and errors) of selected RR Lyrae stars 
are provided in the electronic edition of the Journal; the data format is provided in Tabled! 

In S07 we defined color and variability cuts specifically designed for the selection of 
candidate RR Lyrae stars. We adopt those selection criteria, and expand the probed u — g 
color range from u — g = 0.98 (used in SOT) to u — g = 0.7. This expansion i s motivated by 



theoretical predictions obtained from recent RR Lyrae pulsation models by iMarconi et al. 



(120061 ) ■ They find that, for a given effective temperature, more metal-poo r stars have bluer 



u — g colors. A similar behavior is also detected in main sequence stars (llvezic et al.l 12008 
and references therein) . Thus, the expanded u — g color range might increase completeness 
for more metal-poor RR Lyrae stars; at the same time the sample contamination is not an 
issue because light curves are now well sampled. 

We start the RR Lyrae candidate selection by making the following color 

0.7 <M-^< 1.35 (1) 
-0.15 < ^-r < 0.40 (2) 
-0.15 < r-i < 0.22 (3) 
-0.21 < i-z < 0.25, (4) 

and variability cuts 

• The number of g band observations, Nobs,g > 15 

• The chi-square per degree of freedom in the g and r bands, x^(fl') > 3 and x^(^) > 3 

• The rms scatter in the g and r bands, a{g) > 0.05 mag and cr(r) > 0.05 mag 
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• The ratio of the rms scatter in the g and r bands, 0.7 < a{g)/(j{r) ^ 2.5, and 

• The g band hght curve skewness, 7 ((7) ^ 1, to avoid echpsing binaries (see S07 for 
details). 

Note that the colors in Equations are obtained fro m time-averaged (median) uqr i z mag - 



nitudes, corrected for extinction using the maps from ISchlegel. Finkbeiner fc David (119981 ). 
These criteria are designed to reject stars that are definitely not RR Lyrae, and all but one 
{Nobs,g > 15) should be satisfied by every true RR Lyrae star. The second variability crite- 
rion can introduce incompleteness at the faint end, because photometric errors increase with 
magnitude. We address this incompleteness effect in detail in § 14. 2[ Using these criteria we 
select the initial sample of 3449 candidate RR Lyrae stars (about 1 in 7 are true RR Lyrae 
stars; the fraction of RR Lyrae stars in the full sample is about 1 in 10^), and proceed to 
measure their light curve periods in the next step. 

To estimate the period of pulsati on for candidat e RR Lyrae stars, we use an implemen- 



tation of the Supersmoother routine (jReimanru Il994l ) . The routine uses a running mean or 



running linear regression on the data to fit the observations as a function of phase to a range 



of pe riods, constrained to be between 0.2 and 1.0 days (see Figure 9 in lEyer fc Mowlavi 



20081 ). The main advantage of this routine is that the light curve shape need not be known 
in advance. We supply the Modified Julian date (MJD) of observations, (/-band magnitudes 
and their errors, and the routine outputs the first 5 periods with the lowest chi-square val- 
ues. The ^f-band magnitudes are used instead of u- or r-band magnitudes because they have 
higher signal-to-noise ratios than the tt-band magnitudes, and beca use RR Lyra.e star s have 



a larger amplitude of variability in the g band than in the r band (jivezic et al.ll2000l . S07). 
The precision of periods is about a few seconds, and was determined by comparing periods 
estimated from g- and r-band data. 

Using each of the five periods supplied by the Supersmoother routine, we period-fold a 
(jf-band light curve and then fit to it one of eight ^-b and template R R Lyrae light curves. The 



y-band template light curves were obtained from iLaydenl (Il998l 6 RR type ab templates) 
and K. Vivas (priv. comm., 2 RR type c templates), and are used due to the lack of (7-band 
templates in the literature. The fitting is performed in the least square sense, with the epoch 
of the maximum brightness 0q, peak-to-peak amplitude Ag, and mean magnitude go as free 
parameters. The precision of the epoch of maximum brightness is on the order of a few 
minutes, and was determined by comparing epochs of maximum brightness estimated from 
g- and r-band light curves fitted with l^-band templates. 

To robustly quantify the quality of a template fit we define a x^-like parameter 

( = raedian{\m\g^j.^^^ — mtemplatel/ (^observed)^ (5) 
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where niobserved and aobserved are the observed magnitude and its uncertainty, rritempiate is 
the magnitude predicted by the template, and ^ = 1, Nohs, where N^hs is the number of 
observations. Here we use the median to minimize the bias in C, due to poor observations 
(outhers). The template with the lowest C value is selected as the best fit, and the best-fit 
parameters are stored. 

Initially, we wanted to use C, to separate RR Lyrae stars from non-RR Lyrae stars. Low 
C values, indicating good agreement between the observed light curves (data) and template 
light curves (models), would indicate that a candidate is an RR Lyrae star. Likewise, high C, 
values would indicate that a candidate is not an RR Lyrae star. However, visual inspection 
of light curves initially classified as non-RR Lyrae (high C, values) showed that some of the 
light curves are actually quite RR Lyrae-like, but with a lot of scatter around the template- 
predicted values. We suspect that these stars are RR Lyrae stars that undergo changes in 
their period, l i ght curve shape, or amplitude as a function of time (e.g., the Blazhko effect; 



Blazhko 1907 



KolenberelSosI l By putting a hard limit on C, such RR Lyrae stars would 



be rejected, and the completeness of our sample (the fraction of recovered RR Lyrae stars) 
would decrease. To avoid this, we look for additional criteria beyond the light curve shape 
that can be used to distinguish true RR Lyrae stars from non-RR Lyrae stars. For this 
purpose we have created a "training" set of 479 RR Lyrae stars that were selected by visual 
inspection of their light curves, and use it to explore the variability phase space of RR Lyrae 
stars. Hereafter we refer to the training set as the template-accepted RR Lyrae stars. 

Based on the distribution of template-accepted RR Lyrae stars shown in Figure [21 
we define a line that separates candidate RR Lyrae into template-accepted and ambiguous 
(above the line), and template-rejected (2608 candidates below the line). The ambiguous 
candidates (362 objects) did not have at first visually convincing (RR Lyrae-like) light curves, 
but based on their variability properties (rms scatter in the g, r and g — r) some of them 
might be true RR Lyrae stars. 

The question of whether some of the ambiguous candidates are true RR Lyrae stars is 
especially important at the faint end. As shown in Figure [21 the distribution of template- 
accepted RR Lyrae stars has a sharp cut-off around go ~ 20.5, suggesting a limit on the 
spatial distributi on of halo st a rs. A similar cut-off in the distribution of RR Lyrae stars was 



also discussed by llvezic et al.l (l2000l ). If it can be shown that the ambiguous candidates with 
(?o ~ 20.5 have properties (beside their light curves) that are inconsistent with the properties 
of template-accepted RR Lyrae stars, they could be excluded from the sample, and the claim 
that the halo cut-off is real would become much stronger. 

Aiming to define the RR Lyrae variability space, we now consider how the three groups 
of candidates distribute in the ag vs. 0"^ diagram, where the ag and ar are the root-mean- 



- 8 - 



square (rms) scatter in the g and r bands, respectively. This projection of phase space is 
useful as it separates pulsationally variable stars from other variable sources (S07). As shown 
in Figure [3] (ioj)), the majority of template-accepted RR Lyrae stars exist in a narrow region 
centered on the ag = 1.42(7^ (the rms scatter around ag = 1.420"^ is 0.006 mag) line typical 



of RR Lyrae stars (jivezic et al.l l200d : S07), while the majority of template-rejected and 
ambiguous candidates follow the ag = Or relation that suggests non-pulsational variability. 
It is interesting to note that 13 template-accepted RR Lyrae stars have lower observed Gg 
values (by > 0.03 mag) than what is predicted by the Og = 1.42(7^ (Sa) relation. The g- 
and r-band light curves of two such stars are shown in Figure H] {top and middle). Based 
on the Figure [3] {top) plot, we conclude that the majority of RR Lyrae stars (97%) have 
\ag — lA2ar\ < 0.03, and select all candidates that pass this condition for further analysis. 
The 13 template-accepted RR Lyrae stars found outside the \ag — 1.42(7^1 < 0.03 region are 
also included as we believe they are RR Lyrae stars (based on their light curves), but with 
peculiar variability properties (slightly lower than expected rms scatter in the g band). 

The candidates that pass the \ag — lA2ar\ < 0.03 condition also separate well in the 
ag-r vs. 1^(7^ + 0"^ plot (Figure [3] [middle]), which is sensitive to the covariance between 
the g and r bands. Using the distribution of template-accepted candidates as an indicator 
of where RR Lyrae stars exist in this plot, we make another selection and again plot the 
selected candidates in the g band amplitude vs. mean magnitude plot (Figure [3] [bottom]). 
There are eight template-accepted RR Lyrae stars that do not pass the ag^r vs. ^Jag + 
selection, one of which is shown in Figure H] {bottom). As in the previous paragraph, we 
choose to keep these candidates as they might be RR Lyrae stars with peculiar variability 
properties (higher than expected rms scatter in the g — r color). 

The utility of the (7g_r selection is demonstrated in Figure [5l where we compare g and 
g — r light curves of a type c RR Lyrae (hereafter RRc) and a candidate echpsing binary 
star. Even though their g-hand light curves are quite similar and their periods match those 
of RRc stars, their g — r light curves are quite different. The g — r light curve of the eclipsing 
binary varies much less that that of the RRc star, as the binary is most probably comprised 
of two similar effective temperature stars. On the other hand, the effective temperature of 
the RRc star changes due to pulsations, and the variations in the g — r light curve are much 
greater. 

A comparison of Figures [2] and [3] {bottom) shows that the majority of ambiguous and 
template-rejected candidates do not have variability properties {ag, ar, ag^r) that are consis- 
tent with variability properties of template-accepted RR Lyrae stars. Only 18% of ambiguous 
and 1% of temp late- rejected candidates pass our ag — ar — ag^r variability cuts, while only 
a small fraction of template-accepted candidates (4%) is rejected. Most importantly, al- 



- 9 - 



most all of the ambiguous candidates with go ^ 20.5 were rejected based on their variability 
properties. 

These results suggest that the candidates initially classified as template-accepted are 
indeed true RR Lyrae stars, but that some ambiguous and template-rejected candidates 
might also be true RR Lyrae stars. We have re-examined light curves of ambiguous and 
template- rejected candidates that pass the ag — ar — crg~r variability cuts, and re-classified 
four as template- accepted. The 21 template- accepted candidates that were shghtly outside 
of our selection criteria were also kept, as the visual inspection of their g and g — r light 
curves strongly suggests they are indeed RR Lyrae stars. Finally, we report 483 stars in 
our SDSS stripe 82 RR Lyrae sample, and estimate it to be essentially complete and free of 
non-RR Lyrae stars. 

Our high degree of completeness is achieved by relaxing color and variability cuts from 
S07, and by accepting stars that are slightly outside the selection region if their light curves 
are strongly RR Lyrae-like. We find that the median g — r, r — i, and i — z colors of stars in 
the final sample are well within the ranges defined by Equations and 20 of the selected 
stars (~ 4%) have median u — g color in the 0.85 < u — g < 0.98 range. This finding validates 
the expansion of the u — g range blueward to 0.7 mag (from u — g = 0.98 used in SOT). The 
visual inspection of g and g — r light curves ensures that non-RR Lyrae stars are not included 
in the final sample. 

When time-averaged colors are not available (as is the case for the majority of SDSS 
stars), a range of single-epoch colors might be more useful for the selection of RR Lyrae 
stars. We find that RR Lyrae stars have single-epoch colors in these ranges: 

0.75 < M-^ < 1.45 (6) 

-0.25 <^-r < 0.4 (7) 

-0.2<r-i<0.2 (8) 

-0.3 < z - z < 0.3. (9) 



3. The ugriz RR Lyrae Templates 



The availability of a large sample of RR Lyrae stars with densely sampled light curves 
presents us with a unique opportunity to derive template RR Lyrae light curves for the SDSS 
photometric system. While t emplate RR Lyrae l ight curves for the SDSS ugriz bands do 



exist in the literature (e.g., see lMarconi et al.ll2006l ). they are based on pulsation models, and 



not on actual observations. Recognizing the importance of such templates for future large sky 
surveys that will use the SDSS filter set, such as the Large Synoptic Survey Telescope (LSST; 
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Ivezic et al.ll2008l ) and Pan-STARRS (IKaiser et al.ll2002l ). we now describe the construction 
of template ugriz RR Lyrae light curves. 



3.1. Construction of ugriz RR Lyrae Templates 

Our goal is to construct the smallest set of template ugriz RR Lyrae light curves that 
will model the light curves present in our sample at a level of uncertainty set by photometric 
errors (~ 0.02 mag). To achieve this goal, we use an iterative approach to refine the initial 
set of prototype templates. 

First, for each ugriz band we find RR Lyrae stars that have at least 50 observations 
in a band with photometric errors less than 0.05 mag. The light curves of such RR Lyrae 
are period-folded, and a B-spline is interpolated through the data if the largest gap in the 
phase coverage is smaller than 0.05. The RR Lyrae stars with phase coverage gaps greater 
than 0.05 are not used in the construction of the RR Lyrae templates. The light curves with 
interpolated B-splines (~100 per band) are hereafter referred to as prototype templates, and 
are normalized to have values in the [0, 1] range and maxima at phase = 0. 

In the next step, we fit the prototype templates to each of the 483 light curves and qualify 
a fit as acceptable if C < 1.3. Here it is possible to have more than one acceptable template 
fit per hght curve. The RR Lyrae stars with acceptable fits are grouped by prototype 
templates, and their observed light curves are period-folded and normalized using best-fit 
period, amplitude, mean magnitude, and epoch of the maximum brightness. The normalized 
light curves are grouped on a per template basis, and a B-spline is interpolated through the 
grouped data. With this procedure we smooth out prototype templates by averaging over 
similarly shaped light curves. In other words, the main purpose of the prototype templates 
is to select self-similar subsets of observed light curves. 

Since the number of smoothed prototype templates is about 100 per band, it is possible 
that some of them are redundant and that they may be replaced by a single template without 
decreasing the generahty of the template set. To establish which templates are redundant, 
we fit smoothed templates to RR Lyrae light curves and accept template fits with < 1.3. 
Again, it is possible to have more than one acceptable template fit per light curve. 

We rank the smoothed templates by the number of light curves they fit, and sort them 
in descending order (the highest ranking template fits the greatest number of light curves). 
Starting with the second-ranked template, we count the number of light curves that were 
not fitted by the highest ranking template. We move down the template list and for each 
template count only light curves that were not already fitted by all higher ranking templates. 
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We stop moving down the template list when we reach a template that fits only one light 
curve. 

The templates that survive this pruning are compared with each other, and if there are 
two with maximum difference less than 0.02 (the average photometric error), we keep the 
one that fits a greater number of light curves. Our final, pruned template set contains 12 
u-band (1 is type c), 23 g-hand (2 are type c), 22 r-band (2 are type c), 22 z-band (2 are 
type c), and 19 2;-band templates (1 is type c). The templates are provided in the electronic 
edition of the Journal. The light curve parameters of stripe 82 RR Lyrae stars obtained from 
best-fit ugriz templates are listed Table [2l 

As illustrated in Figure El the set of ~ 20 ab templates per band appears consistent with 
a continuous distribution of light curve shapes. This essentially one-dimensional template 
family can be conveniently parametrized by the phase at which the normalized template 
(with unit amplitude) crosses the value of 0.8 on the rising part of light curve (hereafter 
00.8)- This parameter measures the rise time, as a fraction of the period, and is related to 
other parametrizations of the light curve shape, such as the coefficients of a Fourier expansion 
(see §11221). 

On the other hand, c-type RR Lyrae stars exhibit a strong bimodal distribution of light 
curve shapes in the gri bands, suggesting that this class may include two distinct subclasses, 
with each subclass having very small scatter of light curve shapes. This behavior is further 
discussed in § 14. 6[ 

We also note that our template set is not entirely complete, as illustrated in Figure [71 
There are 15 examples of stars (3%) for which the current template does not describe their 
(unique) light curves at the level of the photometric uncertainty (0.02 mag). 



3.2. A Comparison With Theoretical Template Light Curves 



In FigurelHl we compare empirical RRab templates to theoretical templates from lMarconi et al 



( 120061 ). As the band-pass moves to longer wavelengths (from the g to z band), the light curve 
shape of both empirical and theoretical templates becomes less convex (less curved upward) 
on the descending brightness branch (phase< 0.7). The light curve shape does not change 
significantly between u and g bands in neither empirical or theoretical templates. 

On the other hand, the theoretical templates have a peak after the minimum brightness 
that is not observed in our templates. The differences between our empirical templates and 
theoretical templates are much larger than the errors in our templates. The new pulsation 



models will undoubtedly benefit from constraints provided by the observational template set 
derived here . Additional constra ints may also be obtained from comparisons of theoretical 
(Figure 7 in lMarconi et al.ll2006l ) and observed color-color loops (Figure [HI in this work). 



4. Outer Halo Substructure As Traced by Stripe 82 RR Lyrae Stars 

In this section the spatial distribution of RR Lyrae stars is used to trace the substruc- 
tures in the outer halo. First, we estimate the heliocentric distances of RR Lyrae stars using 
their best-fit g- and r-band light curves, then estimate the sample completeness as a function 
of magnitude in Section 14. 2[ The observed spatial distribution of RR Lyrae stars is used 
to compute their number density distribution, which is then compared to a model number 
density distribution. This comparison reveals an abrupt change in the number density dis- 
tribution at a distance of ~ 30 kpc, suggesting a transition from the inner to outer halo, as 
well as the presence of a prominent halo overdensity at 80 kpc. 



4.1. Distance Estimates 



The heliocentric distances, d, for RR Lyrae stars are estimated from measurements of 
their mean dereddened V magnitudes, (V^), and by assuming an average value for their 
absolute magnitudes. My 

d = 10((V'>-Mv.+5)/5 (IQ) 

For the absolute magnitude in the V band we adopt My = 0.6, calculated using the 
Cacciari k Clementinil J2003h My - [Fe/H] relation 



Mv = (0.23 ± OM)[Fe/H] + (0.93 ± 0.12) 



'11^ 



and adopting [Fe/H] = —1.5 as the mean halo metallicity (jivezic et al.ll2008l ) 



The values of {V) are determined by integrating best-fit template y-band light curves 
after transformation to flux units. The best-fit template V"-band light curves are synthe- 
sized from the best-fit g- and r-band tem plates (qr,r[ (l )) and rLc{(p)) using a photometric 
transformation given by Equation 10 from 



Ivezic et all (l2007f ): 



9rLc{4>) = 9lc{4>) - rLc{4>) (12) 
V^Lc(0) = 9Lc{(t>) + O.O688(?rLc(0)' - O.2O56(?ric(0)' - O.3838(7ric(0) - 0.053 (13) 



where G [0, 1) is the phase of pulsation, and where gici.'P) and rici.'P) are corrected 
for interstellar extinction. We have also established that the measured light curves are 
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aligned in phase among the bands to within the measurement errors (to within 1% in phase). 
These "flux-averaged" magnitudes are used because pulsation models show that they better 
reproduce the behavior of "static" magnitud es (the magnitudes the stars would have if 
they were not pulsating, iMarconi et al.ll2006l ). The uncertainty in (y) is mainly due to 
photometric errors in the g and r bands, and is estimated to be ~ 0.03 mag. 



Following the discussion by lVivas fc ZinnI (120061 ) (see their Section 4), we adopt cjm, 



'^^M^^^^Y + {^MyY = 0.1 mag as the uncertainty in absolute magnitude, wh ere (t^Mv^^ 



0.07 mag is the uncertainty due to metallicity dispersion {(T[Fe/H] = 0.3 dex, llvezic et al. 
20081 ). and /S.My = 0.08 mag is the un certainty due to RR Lyrae evolution off the zero-age 
horizontal branch (jVivas fc ZinnI |2006| ). The estimated fractional error in the heliocentric 
distance [d/a d) of 0.06 is smaller, by a factor of ~ 2 — 4, than the fracti onal distance error 
for K giants f|Dohm- Palmer et al.l 120011) . M gi ants (IMajewski et al.l 120031 ) or main sequence 
stars ( Newberg et al. 20021 : Juric et al. 2008). The hehocentric distances and equatorial 
coordinates of stripe 82 RR Lyrae stars are listed in Table [3l 



4.2. Sample Completeness as a Function of Magnitude 

The fraction of recovered RR Lyrae stars (the completeness) as a function of apparent 
magnitude is an important question that needs to be considered before the spatial distri- 
bution of RR Lyrae stars can be discussed. We generate 1000 RRab synthetic gf-band light 
curves with 30 data points, amplitudes of 1.0 mag, periods of 0.6 days, and with flux- 
averaged magnitudes of {g) = 21.6 (not corrected for ISM extinction). An RR Lyrae star 
with {g) = 21.6 would be located at ~ 120 kpc, assuming Mg = 0.6 and 0.6 mag of ISM 
extinction in the g band (the highest observed extinction in stripe 82). Typical random pho- 
tometric errors were added as a function of magnitude using Figure 1 from SOT. The periods 
for the synthetic light curves were estimated using Supersmoother, and the amplitudes and 
mean magnitudes were estimated from best-fit templates. 

We find that the periods are recovered with 0.05% (or better) accuracy, and flux- 
averaged magnitudes with 0.03 mag (or better) accuracy. Similar results are obtained for 
RRc stars (the amplitude of synthetic light curve is 0.5 mag in that case). If we define an 
RR Lyrae as successfully recovered if its fiux-averaged magnitude is estimated with 0.03 mag 
accuracy, then the stripe 82 RR Lyrae sample is essentially 100% complete to 120 kpc, even 
in regions where the extinction is as high as 0.6 mag {R.A. > 50° and R.A. < 308°). 

Even though both RRab and RRc types have the same completeness as a function of 
apparent magnitude, we only use RRab types when studying the spatial distribution of RR 
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Lyrae stars, because they represent a homogeneous and robust sample with a well-studied 
absolute magnitude vs. metallicity distribution. A repeat of the analysis presented below, 
but with RRc stars included, did not change our conclusions based exclusively on RRab 
stars. 



4.3. Observed Spatial Distribution 



Since SDSS stripe 82 is quite narrow in declination (2.5° compared to 118° in R.A.), we 
assume Dec = for all stars, and describe their positions in the R. A.— heliocentric distance 
space. The Cartesian coordinates of stars in this space are defined as 



X = dcos{a) 
y = dsm{a), 

where d is the heliocentric distance in kpc, and a is the right ascension angle. 



(14) 
(15) 



The spatial distribution of 366 RRab stars in the R. A.— heliocentric distance space 
is shown in Figure dUl For comparison, we also show the distribution of RR Lyrae stars 
randomly drawn from a smooth mode l number de r isity d istribution, based on an analysis of 
main sequence stars within 20 kpc by lJuric et al.l (120081 ) 

Rq 



Pmodel 

{R,Z) 



^R^ + {Z/quf 



(16) 



where p^^ 



4.2 kpc is the number density of RRab stars at Rq = 8.0 kpc (IVivas fc Zinn 



20061 ). qn = 0.64 indicates the halo i s oblate ( "squash ed" in the same sense as the disk), and 
nn = 2.77 is the power-law index (jJuric et al.l l2008l ). The expected number of RR Lyrae 



stars (according to p^^ei) ^^e 5 to 120 kpc range of st ripe 82 is 302. The value of 
within ±0.5 kpc~^ from va 



^RR ; 



Suntzeff. Kinman fc Kraft 



ues reported by other authors (jPreston. Shectman &: Beerslll991 



fil- 



ls 



199ll : IWetterer &: McGrawlll996l ). and the uncertainty in uh and 



<iH is < 0.2 (jJuric et al.ll2008l ). We also repeated the following analysis using a spherical. 



1/R^ halo model and obtained similar results and conclusions (i.e. detection of significant 
overdensities and the steepening of the halo density profile). 

The two distributions shown in Figure [10] appear remarkably different. The observed 
distribution, shown in the top panel, has two overdense regions within 30 kpc; the one at 
R.A. ~ 30° c an be identified as the trailing arm of the Sagittarius dwarf spheroidal (dSph) 
tidal stream (jivezic et al.l 120031) . while the othe r at i?.y4 ~ 330° is likely associated with the 



Hercules-Aquila Cloud flBelokurov et al.l l2007bl ) . Beyond 30 kpc, the observed distribution 



is quite inhomogeneous in the R.A. direction, with the majority of RR Lyrae in the 330° to 
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0° region, and only a few RR Lyrae stars outside of it. On the other hand, the distribution 
of stars generated using Equation [16] looks more homogeneous, as expected. 



4.4. Bayesian Number Density Maps 

We com pute the loca l num ber density of RR Lyrae stars using a Bayesian estimator 



developed by llvezic et al.l (l2005l . see their Appendix B2 and Equation B6). The Bayesian 



metric presented there can be recast in a much simpler form 



^stars y^,j-l 



where Nnn is the number of nearest neighbors to which the distance r in a /c-dimensional 
space is calculated, and Ngtars is the number of stars. The dimensionless constant C involves 
complex integrals, but can be easily determined empirically once the geometry of the space 
in which densities need to be calculated is known. 

To estimate the constant C we use the following procedure. We populate the stripe 
82 region with Ngtars = 302 stars randomly drawn from the smooth model number density 
distribution (EquationHH]), and calculate the Bayesian metric on a 0.5x0.5 x—y grid assuming 
Nnn = 8, k = 2, and C = 1. The nearest neighbors are fou nd using an implementation of 



the nearest neighbor searching algorithm (lArya et al.l 119981 ). The number density is then 
calculated as 

pix,y) = pBayesix,y)/d{x,y), (18) 

where the Bayesian metric is divided by the heliocentric distance to account for the increase 
in volume of a grid point at greater heliocentric distances. We generate 1000 p(x, y) density 
maps, and average them out to obtain pavgix, y) = ^ p{x, ?/)/1000. The constant C = 1030 is 
estimated by requiring the median value of Pavg{x, y) to equal the median value of Pmodeii^y u)- 

The sparseness of our sample (366 RR Lyrae in a volume of 38,000 kpc^) will result 
in random (Poisson) number density fluctuations even for a smooth underlying density dis- 
tribution. To quantify the scatter in number density maps and to estimate the fraction of 
highly overdense and underdense pixels (| log(p/p^o^g;)| > 1) due to Poisson fluctuations, we 
generate 1000 random RR Lyrae star samples of the same size as the observed sample by 
drawing them from the smooth model, and compute ^og{p/p^^^i) maps for each sample. We 
find that the median rms scatter in smooth model-based maps is 0.3 dex (or about a factor 
of two on the linear scale), and that the fraction of highly overdense and underdense pixels 
(I log(p/p^^^gJ| > 1) is only 10~^. We conclude that the precision of the Bayesian density 
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estimator, given the sparseness of our sample, is about 0.3 dex, and that highly overdense 
and underdense pixels are not likely to be caused by Poisson fluctuations. 

With the Bayesian density estimator now fully defined, we compute the number density 
y) on a X — y grid for the observed distribution of RR Lyrae stars and for the sample of 

stars randomly drawn from the smooth model distribution. The computed number density 

maps are divided by Pmodeii^^y)'! results are shown in Figure [TTl 

Figure [TT] (top) vividly shows the striking difference between the observed number den- 
sity distribution of RR Lyrae stars and the smooth model number density distribution (Equa- 
tion [16]). While the density fluctuations (log(p/p^^^g;)) in the model-based map (Figure [TT] 
[bottom]) are consistent with Poisson noise (amplitudes are < 0.3 dex), the observed (top) 
map has several overdense and underdense regions (by a factor of 10 or more relative to 
Pmodei)- The Sagittarius dSph tidal stream (trailing arm) is visible at ~ 25 kpc and 30°, the 
Hercules-Aquila cloud is located at < 25 kpc and ~ 330°, and a previously detected, and 
now confirmed, structure can be seen at ~ 80 kpc and ~ 355° (the "J" overdensity in S07 
Figure 13). The fraction of pixels with | ^og{p/ p^^^i)\ > 1 in the observed number density 
map is about 0.3 (dominated by underdense pixels), or 5 orders of magnitude higher than 
expected for Poisson fluctuations. 

Figure [TT] (top) also shows that some overdensities detected by SOT are not reproduced 
in this work. These overdensities were caused by non-RR Lyrae stars in the S07 candidate 
RR Lyrae sample (Table [1] lists the number of RR Lyrae and non-RR Lyrae stars in each 
SOT overdensity). About 70% of S07 RR Lyrae candidates are true RR Lyrae stars, with 
the re-classified non-RR Lyrae stars dominated by 6 Scuti stars (based on their light curves 
and periods). These results are in good agreement with S07 estimates. The overdensities 
labeled D, F, H, I, K, and L now have lower numbers of true RR Lyrae stars compared to 
S07, and are no longer significant overdensities. These overdensities were also speculated in 
S07 to be false detections. The significance of the E, G, and M overdensities decreased as 
the number of true RR Lyrae stars decreased compared to S07, and no longer seem to be 
real substructures. 

The extension of the J overdensity in S07 in the radial direction due to the uncertainty 
in distance (the "finger of God" effect) is now smaller than in S07. In S07, t he apparent 



V^-ban d magnitudes of RR Lyrae stars were computed using Equation 3 from llvezic et al. 



(I2OO5I ) and using single-epoch measurements: 

Vs%t = r - 2M{g - r) + 0.355. (19) 



While this approximation works fairly well for RRab stars (the rms scatter around the fit is 
~ 0.1 mag), it produces an offset of ~ 0.3 mag when applied to RRc stars, due to a 0.1 mag 
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offset in their g — r colors (RRc stars ha ve bluer q — r color s than RRab stars). The higher 
fraction of RRab than RRc stars i n the llvezic et al.l (120051 ) RR Lyrae sample may explain 



the resulting bias in Equation [12] (llvezic et al 



2005 



did not fit Equation [12] separately for 
RRab and RRc stars). In this work we use flux-averaged \^-band apparent magnitudes and 
obtain more precise distance estimates for both RRab and RRc stars. 

The structure at ~ 80 kpc and ~ 355° (labeled J in S07) is most probably a tidal stream, 
although the possibility that it is a dSph satellite galaxy cannot be excluded with the data 
presented here. It is unlikely that this structure is associated with the S agittarius dSph 
stream because it is ~ 22 kpc away from the best-fit orbital plane defined by lMajewski et al. 



(|2003l). For compari son, the Sgr dSph trailing tidal stream and th e stream detec t ed by 



Newberg et al.l (120031 ) are both within 2 kpc of the Sgr dSph plane. IWatkins et al.l (120091 ) 
point out that the location of the J overdensity is close to the Magellanic Plane, though 
at a larger galactocentric distance than the Clouds. They also proposed to name the J 
overdensity as "the Pisces overdensity", following the tradition of naming stellar streams 
and structures according to the nearest constellation. We adopt their proposal hereafter. 

The width of the Pisces overdensity /stream in the radial direction (~ 4 kpc) is consistent 
with the scatter due to uncertainty in distance, while the scatter in the R.A. direction (~ 4 
kpc) might be a projection effect (a combination of the stream's incidence angle and thickness 
of the stripe 82 plane), or evidence of the stream's width in the R.A. direction. To illustrate 
the scale and position of this and other detected structures, we provide an animation showing 
a "fly-by" of the stripe 82 region. The animation is provided in the electronic edition of the 
Journal; a single frame is shown as Figure j 



4.5. The Steepening of the Halo Density Profile 

The observed number density map (Figure [11] [top]) suggests that the smooth model 
(Equation [TB]) systematically overestimates the observed number densities at galactocentric 
distances Rqc > 30 kpc. We investigate this discrepancy in Figure [13] (top) by comparing 
the observed log(p/p^^^e«) distributions within Rgc < 30 and Rgc < HO kpc with the 
distribution obtained from model-based number density maps. After a 33% increase in the 
model normalization {pq^ from Equation [T6]) . the observed and model-based distributions 
match well for Rqc < 30 kpc. About 80% of the halo volume within Rgc < 30 kpc 
is consistent with the smooth model (within Poisson noise limits, \^og{p/ p^^^^i)\ < 0.3); 
about 20% of the halo volume is found in predominantly overdense, substructures. The 
overdense tail of the distribution is more pronounced beyond 30 kpc. At these larger radii, 
the cumulative distribution (Figure [12] [bottom]) indicates that most of the probed halo 
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volume is underdense, relative to the adopted smooth model (which is extrapolated from 
observations of main sequence stars within ~20 kpc). 

The apparent need for a steeper density power-law index at Rgc > 30 kpc is also il- 
lustrated in Figure [TU where we plot the number density of RR Lyrae stars as a function 
of position angle and galac tocentric distance. F or Rgc > 30 kpc and directions that avoid 



the Pisces overdensity, the IWatkins et al.l (120091) best-fit "broken" power-law (based on es- 



sentially the same RR Lyrae sample) provides a good fit with a ri/^ = 4.5 power-law slope 
(apart from normalization: we had to scale down their normalization, taken from their arXiv 
posting, by a factor of 10). 



The apparent good agreement between the data and smooth, lJuric et al.l (120081 ) halo 
model in the 340° < R.A. < 0° direction at ~ 80 kpc (Figure HH [bottom left]) is due to 
Pisces overdensity. 

To quantify the statistical significance of the underdense regions relative to a smooth 
model, we have divided stripe 82 into 6 regions, and for each computed the expected number 
of RR Lyrae stars by integrating p^odei- The statistical significance of the underdense regions 
is calculated by finding the absolute difference between the observed (Nobs) and expected 
(Nexp) number of RR Lyrae stars (based on the oblate halo model. Equation 1T6|) . then 
dividing it by the Poisson uncertainty ( N^xp) ■ The results are presented in Table ISj the 
underdense regions are statistically significant at the (1 — 3)a level compared to the oblate 
halo model (Equation HE]). These underdensities cannot be explained by incompleteness as a 
function of apparent magnitude (since there are dozens of RR Lyrae detected up to 100 kpc 
in the Pisces stream region) or angle (because SDSS provides very homogeneous coverage of 
the sky). The simplest interpretation is that the spatial profile steepens at ~ 30 kpc. This 
hypothesis is also supported by the distribution of halo main sequence stars, discussed in 
Section 



Using the selection criteria from Section 14.61 we have divided the sample of RRab stars 
into Oosterhof I (Oo I, 289 stars) and Oosterhof II (Oo II, 77 stars) subsamples and compared 
their radial profiles in Figure ITS) (top). In the direction that does not overlap with the Sgr 
tidal stream and the Hercules-Aquila Cloud (340° < R.A. < 30°), and for Rqc < 30 kpc, 
the Oo I subsample shows a shalower profile than the Oo II subsample. The best-fit power 
law indices are —2.3 and — 4.5, respectiv ely. The fit for the Ool subsample is in excellent 



agreement with the result by lMiceli et al.l (l2008l ) obtained over a much larger sky area (1430 



deg^ vs. ~300 deg^). For the OoII subsample, our best fit is steeper (Miceli et al. obtained 
a best-fit index of —2.9). The number of stars with Rgc > 30 kpc is too small (only 23 Oo 
II stars) to robustly measure the profile difference beyond the 30 kpc. We also note that 
the radial number density profiles of RRab and RRc stars have consistent best-fit power law 
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indices (—2.9), as shown in Figure [15] [bottom). Interestingly, the Hercules-Aquila Cloud 
seems to be dominated by Oo I stars; the Pisces stream includes only Oo I stars. 



4.6. Period- Amplitude Distribution 



As suggested by ICatelanI (120091 ). the period-amplitude distribution of RR Lyrae stars 



may hold clues about the formation history of the Ga lactic halo. Catelan points to a sharp 



division (a dichotomy first noted by lOosterhofI Il939l ) between Oosterhof type I (hereafter 



Oo I, average periods of RRab stars (Pab) ~ 0.55 days), and Oosterhof type II (hereafter 
Oo II, average periods of RRab stars {Pat) ~ 0.65 days) globular clusters, with only very 
few clusters in between. On the other hand, the dSph satellite galaxies and their globular 
clusters fall preferentially on the "Oosterhof gap" (0.58 < {Pab) < 0.62, see his Figure 5). 
In light of these findings, it is interesting to examine the period-amplitude distribution of 
stripe 82 RR Lyrae stars in general, and the position of RR Lyrae stars from the Sgr and 
Pisces streams in particular, in the period-amplitude diagram. 

The period-amplitude (or Bailey) diagram for RR Lyrae stars listed in Table [2] is shown 
in Figure [161 The stars are separated by the "quality" of their (/-band light curves. A star 
with a low quality light curve (a lot of scatter around the best-fit template, ( > 2.3) may 
indicate a Blazhko variable RR Lyrae star. The data points are color-coded by the median 
g — r color at minimum brightness {{g — r)min) and by a template parameter 0Qg (which 
measures the rise time). 

We find that {g — r)mm of RRab and RRc stars becomes redder (effective temperature 
at minimum brightness decreases) as the period of pulsation increases (cooler stars pulsate 
more slowly). The (7-band light curve shape (i.e., the best-fit template) for RRab stars 
correlates strongly with the amplitude (and period), becoming on average more asymmetric 
(0Qg increases) with increasing amplitude (and decreasing period). Using a linear fit, we find 
that the expression 

0^.8 = 0.961 + ^m^Ag - 0.249P (20) 

reproduces the template measured with an rms scatter of 0.025 (the measurement error 
for 0Qg is ~0.01). The distribution of light curve shapes for RRab stars appears to be 
continuous. 

The shape of RRc stars, on the other hand, appears to be bimodal, a finding supported 
by the light curves shown in Figure [T71 We emphasize that this conclusion is independent of 
template fitting; the only information needed to fold the light curves in this figure is the best- 
fit period. The RRc stars with more asymmetric shapes ((^gg ~ 0.8, dashed line template 
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in Figure [6] [bottom]) also have on average higher gf-band amphtudes than RRc stars with 
more symmetric curves (by ~ 0.1 mag). These trends for RRab stars have been known fo r 



some time (in the Johnson UBVRI photometric system; ISandage. Katem &: SandagdllQSlI ). 
while the strong bimodality in the (7-band light curve shape of RRc stars is apparently a new 
result. 

The period-amplitude distribution of RRab stars in Figure [16] appears to be bimodal, 
with the two groups possibly being Oosterhof type I and type II stars. We fit a quadratic 
line 

Ag = -3.18 - 26.53 log{Period) - 37.88 \og{Periodf (21) 
to hig h quality light curve stars in the ma in locus, and note that the fitted line is very similar 



to the ICacciari. Corwin fc Carney! (120051 ) period-amplitude line of M3 globular cluster RR 
Lyrae stars. Since M3 is the prototype Oo I globular cluster, we label the main locus as 
the Oo I locus. The stars with low quality light curves are not used in this fit, as they are 
probably Blazhko variables that might not be observed at their maximum amplitude, and 
therefore tend to scatter towards lower amplitudes in the Bailey diagram. We calculate the 
period shift, Alog{Period) (the log{Period) distance at a constant amplitude from the Oo 
I locus line) as 

Alog{Period) = log{Period) — log{Period interpolated from the Oo I locus) . (22) 

The distribution of A log{Period) values is centered on A log{Period) = by definition (the 
position of the Oo I locus), and has a long-period tail. Even though we do not see the clearly 
di splaced secondary peak that is usually associated with an Oo II component (see Figure 21 



m 



Miceli et al.ll2008l ). hereafter we will refer to stars in the long-period tail as Oo II stars. 
Assuming that the distribution around A\og{Period) = is symmetric for Oo I stars, we 
define RRab stars with A\og{Period) < 0.03 as Oo I, and those with A log(Perio(i) > 0.03 
as Oo II. 

The period- amplitude distributions of RR Lyrae stars in the Sgr stream (15° < R.A. < 
44°, 18 < d/kpc < 35) and the Pisces stream [R.A. > 345° R.A. < 0°, 75 < d/kpc < 100) 
are shown in Figure [13 The Sgr stream has nine RRab Oo II stars while the Pisces stream 
has none; if the Pisces stream had the same fraction of Oo II stars as the Sgr stream, four 
would be expected. We have already discussed the differences in the spatial profiles between 
Oo I and Oo II stars in Section 14.51 and provide evidence that Oo II stars have metallicities 
about 0.3 dex lower than Oo I stars in the next Section. 
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4.7. Metallicity Estimates for RR Lyrae Stars 



In order to obtain a rougli estimate of tlie metallicity of the Pisces and Sagittarius 
stream RR Lyrae stars, we consider several methods. The most accurat e metallicity esti- 
mates for RR Lyrae stars (±0. 2 dex, Vivas et al. 2008 : Prior et al. 20091) ca.n be obtained 
by employing the AS* method (jPrestonlll959l ). or one of its variants (lLaydenlll994l ). on the 
observed RR Lyrae spectra. In the absence of spectra, sever al photometr ic metallicity in- 
dices can be used, such as the ultr aviolet blanketing index (ISturchI Il966l). the Stromgren 



mi index ( Epstein fc Epstein 1973 ). or the Walraven system (jBlancol 



1992 and references 



therein). These photometric metallicity indices are less accurate than the spectroscopic es- 
timates (±0.5 dex vs. ±0.2 dex), and are usually used when the stars of interest are faint, 
or when samples of RR Lyrae stars are large. Another method to estimate metallicity for 
RR Lyrae stars is based on the light curve shape, as discussed below. 



4.7.1. Photometric Metallicity 



We now attempt to develop a photometric metallicity index similar to Sturch's ultravi- 
olet blanketing index, but based on the SDSS ugr bandpasses (instead of the Johnson UBV 
bandpasses). The basic idea behind this index is that the correlation between the U — B 
and B — V color at minimum brightness can be used to estimate the metallicity of a RRab 
star. The physics of this estimate is simila r to that of photometric metallicity method used 
for main sequence stars (jivezic et al.ll2008l and references therein). 



In Figure [T9l (top), we investigate if such a correlation exists between the u — g and g — r 
colors on a sample of 51 RRab stars with measured spectroscopic metallicities. The metal- 
licities were obtained from SDSS spectra with the pulsation phase taken into account, and 
were pr o cesse d through the SEGUE Spectral Parameters Pipeline (SSPP; iLee et al.ll2008l ) by 
De Led ( l2008l ). separately from the SDSS Data Release 7 spectra. The random uncertainty 
in the De Lee et al. spectroscopic metallicity is estimated at ~ 0.3 dex, with systematic 
uncertainties of a similar magnitude. We have compared the De Lee et al. metallicity values 
to phase-averaged values listed as part of the SDSS Data Release 7 and found a median 
offset of —0.15 dex and rms scatter of 0.2 dex. 



The conclusion derived from data shown in Figure [19] is that the most metal-poor 
RR Lyrae {[Fe/H] < —2.0) usually have {u — g)min < 1.11- These metal-poor stars are 
also predominantly classified as Oo II stars, as shown in Figure [T^ {bottom). The median 
{u — g)min of RRab stars in the Pisces and Sagittarius streams is ~ 1.13 mag, resulting in a 
weak constraint that their metallicities are in the range —2.0 < [Fe/H] < —1.5] (on the De 
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Lee metallicity scale). 

Motivated by the [Fe/H] vs. u — g color distribution shown in the bottom panel of 
Figure Uni we tried to estimate the photometric metallicity using only the {u — g)min color. 
We first compute the median spectroscopic metallicity in 0.02 wide steps of the {u — g)min 
color, from {u — g)min = 1-07 to {u — g)min = 1.25, then fit a parabola to these medians (see 
Figure [20]) . The best-fit expression is 

[Fe/HUoto = -46.47 + 72.36(^i - g) - 28.98(m - g)l,^. (23) 

In the validity range 1.05 < {u — g)min < 1-20, the best-fit values vary from —2.4 to —1.4; 
the largest value of the derivative of the best-fit is 0.1 dex/0.01 mag for {u — g)min = 1-07. 
Of course, the absolute metallicity error (systematic error) in this expression is inherited 
from the input data provided by De Lee et al.; it could be as large as several times 0.1 
dex. We emphasize that this estimator should not be applied outside the color range 1.05 < 
{u — g)min < 1-20, and that the large scatter of the calibration sample around the best-fit 
relation (~ 0.3 dex) should be noted. 

Equation [23] predicts [Fe/H]photo ~ —1.7 de x for the metallicity of both the Sgr and 



Pisces streams. A similar value was obtained by IVivas. Zinn fc GallartI (120051 ) from spec- 
troscopic observations of RR Lyrae stars in the leading arm of the Sagittarius tidal stream. 
However, the Equation [2S1 result for the trailing arm of the Sgr stream (discussed here) 
is lower by 0.5 dex than its metallicity obtained using halo main sequence stars (in Sec- 
tion [5lT] below). A possible explanation for this discrepancy may be a systematic offset in 
the SSPP metallicity scale for hotter stars. For 45 RR Lyrae from the Sgr tidal stream 
(15° < R.A. < 44°, 18 < d < 35 kpc), we find using the SDSS DR7 SSPP data a median 
radial velocity of —150 km s^^ (with an rms scatter of 35 km s~^, indicating a fairly cold 
stream), and a median metallicity of [Fe/H] = —1.5 (rms scatter of 0.2 dex). If the true 
median metallicity of the Sgr tidal stream (trailing arm) is [Fe/H] = —1.2, as suggested in 
Section \5A] then the SSPP underestimates the metallicities of RR Lyrae (and possibly blue 
horizontal-branch) stars by 0.3 dex. The combination of this —0.3 dex offset and a —0.15 
dex offset between the DR7 and De Lee metallicities could explain the 0.5 dex difference 
between the median metallicity based on Equation [23] and the main sequence-based value 
obtained in the next section. The SSPP results for hotter stars are presently being refined, 
based on comparisons with high-resolution spectroscopy, and this may resolve the apparent 
discrepancy in the near future. 
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4-7.2. Relationship between Metallicity and Light Curve Shape 

Another potentially useful method for estimatin g the metallicity of RR Lyrae stars is 
based on the Fourier decomposition of light curves (jJurcsik fc Kovacs 119961) . Their resu lt 



that the Fourier parameters are correlated with the metallicity (see also lKunder fc Chaboyer 
20081 for an independent confirmation of this correlation) implies that our templates should 
also correlate with metallicity. Furthermore, they have established that the expansion coef- 
ficients are well correlated among themselves, which is consistent with our finding that light 
curves are essentially a one-dimensional family. 



Following ISimon fc ClementI (119931 ). I Jurcsik fc Kovacd (119961 ) parametrized their metallicity- 
light curve shape relation using the phase of the third harmonic (for the sine expansion), 
031. We have expanded our type-afe templates using six sine terms (the fitting residuals for 
five terms are still comparable to our template errors, and more than six terms does not 
result in substantial best-fit improvement) and found a nearly perfect relation between our 
(7-band template light curve parametrization (^gg) and the Fourier expansion (03i): 



()3i = -6.640^.8 + 10.81 



(24) 



which reproduces directly measured 03i (from the Fourier expansion) with a rms of 0.07 (a 
typical fractional error of 1-2%), roughly the same as the error in 03i due to the template 
error for (p^^ (~0.01). Examples of the best Fourier fits to our templates are shown in 
Figure [211 It is evident that the third harmonic controls the position of the light curve 
minimum, which in turn determines the value of ^gs- 



With this relation, we ca n now compute the metalhcity estimates using the ISandage 



(I2OO4J ) expression (given on the lZinn fc Westlll984j scale: the Jurcsik & Kovacs scale is about 
0.3 dex more metal rich): 



[Fe/H]sanOA = -5.49 - 5.664P + 1.412031 



(25) 



This expression is valid for P > 0.43 days and 4.7 < 03i < 5.7, and has a precision of 0.2 
dex. The difference between the [Fe/HjsanOi and the SDSS spectroscopic metallicity for 51 
stars with the latter data has a median offset of -0.07 dex and rms scatter of 0.33 dex. Given 
that the errors in SDSS values are about 0.3 dex, the agreement is satisfactory. Indeed, it is 
encouraging that the median offset is so small given that the methods and their calibration 
are entirely independent. 

Using the SDSS metallicity for 51 stars, we have also determined a best-fit linear com- 
bination of period, g-haiad amplitude and 03i 



[Fe/H] 



SDSS 



-0.66 - 3.65P - O.mAg + 0.313031 



(26) 
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which reproduces the spectroscopic metallicity with an rms scatter of 0.26 dex. The scatter 
is smaller than for [Fe/ H]sanO'i because Equation [26] is directly fit to the data. On average, 
[Fe/ H]sDss is 0.1 dex higher than [Fe/ H]sanQ4:^ with an rms scatter of their difference of 
only 0.13 dex; hence, the two expressions are practically equivalent. Equation [221 together 
with Equation [2ll implies that for fixed values of period and amplitude, the phase of the 
light curve minimum increases as the metallicity decreases; on average, low-metallicity RR 
Lyrae have shorter rise time than metal rich RR Lyrae. 

With the light curve based metallicity given by Equation [261 we revisited the [Fe/H] 
vs. u — g relationship given by Equation [23] using a six times larger and much less selection- 
biased sample. The results show a much smaller gradient of the u — g color with metallicity, 
about 0.02 mag/dex, and call in question the reliability of the photometric metallicity esti- 
mator given by Equation [23] 

Using metallicity estimates computed by Equations [23] and [25] we have compared the 
metallicity distribution of RR Lyrae in the Sgr and Pisces streams to the rest of the sample. 
To within 0.05 dex, there is no detectable difference. On the other hand, the light curve 
based metallicity distributions for Oo I and Oo II stars are different. The median metallicity 
for Oo I stars is 0.31 dex higher when using Equation [23] and 0.46 dex higher when using 
Equation [25] However, both the separation of Oo I from Oo II stars, and the light curve 
based metallicity use the information about period and amplitude (e.g., a large period results 
in Oo II classification and low metallicity). A somewhat smaller, but still significant (~ 2cr) 
difference in the median metallicity of 0.2 dex is found between 37 Oo I stars and 14 Oo II 
stars with available spectroscopic metallicities. 

We have also tested hght curve-based metallicity estimates using a sample of 11 RR 
Lyrae s tars with accurate light curves, and for which there are reliable spectroscopic metal- 



licities ( IVivas et al.l 120081 ). We fit our ^f-band template set to their light curves, and apply 
Equations [25] and [26] to compute [Fe/H]. The differences between these light curve-based 
metallicity estimates and spectroscopic values show an rms scatter of only 0.14 dex. The 
median offset is 0.23 dex, with the spectroscopic values being higher. This offset may in- 
dicate a need to recalibrate the Sandage relation for use with our templates, because both 



his rel ation and the spectroscopic metallicities are supposed to be tied to the IZinn &: West 



(119841 ) scale. 



In summary, given the strong correlations between independe nt light curve quantities 



peri od, amplitude, and the shap e parameter 03i, and the results of [Jurcsik fc Kovacd (119961 ) 
and lKunder fc Chaboyerl ( 20081 ). it is likely that RR Lyrae light curves contain information 
about metallicity. However, the currently available spectroscopic subsample for the SDSS 
stripe 82 region is insufficient to convincingly demonstrate this relationship. 
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We note that IWatkins et al.l (120091 ) . using essentially the same data set as we did, have 
claimed that the metallicity for RRab stars can be determined from SDSS light curves with 
an rms of 0.25 dex by applying the Jurcsik & Kovacs method, in apparent conflict with our 
conclusions. A possible explanation for this discrepancy is that they fit a large number of 
free parameters (eig ht) to small sarnples o f stars (~100). We have used the same functional 
form as proposed by lWatkins et al.l (120091 ) to fit ran^iom/?/ generated values for [Fe/H] and 
other parameters. In 20-30% of simulated cases, we o btain very similar dist ributions in the 
[Fe/H] Jit vs. [Fe/ H]input diagrams as those shown by lWatkins et al.l (120091 ). 



5. The Distribution of main sequence Stars in Stripe 82 

The strength of the evidence for a steepening of the halo stellar number density profile 
discussed in § 14. 51 is hmited by the small size of the RR Lyrae sample, as well as by difficulties 
associated with using RR Lyrae stars as tracers (see Section [T]). A more robust case can 
be made by tracing the distribution of t he much mo r e num erous main sequence stars, using 



an approach such as that discussed by lJuric et al.l (120081 ). Over most of the SDSS area 



the distribution of blue main sequence stars can be mapped out to a distance limit of ~20 
kpc, which is insufficient to test the transition at ~30 kpc. However, the co-added stripe 82 
data (Annis et al., in prep) have a limiting magnitude about two magnitudes fainter than 
single-epoch SDSS data, and can be used to map the number density distribution of blue 
main sequence stars out to ~40 kpc. 

Figure [22] shows the distribution of stars, in several narrow color ranges, as a function 
of apparent magnitude and position along the stripe 82. In each color bin, stars have 
approximately the same absolute magnitude, thus more distant stars are fainter. The data 
clearly show a well-defined isolated overdensity of stars at R.A.~ 20° — 55°. As the color 
becomes redder the overdensity is detected at fainter magnitudes, ranging from r ~ 21 — 22 
for the bluest bin and disappearing at the faint end (r ~ 23.5) in the reddest bin. This 
behavior is consistent with the behavior of the main sequence stars, as illustrated in the top 
middle panel in Figure! 



The top right and bottom panels in Figure [231 exhibit qualitative evidence for the steep- 
ening of the halo density profile; in the color range 0.3 < g — i < 0.5, and away from the 
Sgr stream, the differential counts of stars significantly decreases (by about a factor of two) 
at magnitudes r > 22. If the density profile followed the canonical power law, the 
differential counts would be approximately constant. In order to quantitatively address this 
effect, we co mpare the observe d counts to a simulated sample based on a smooth oblate halo 



model from jjuric et al.l (120081 ). As shown in Figure [2H the model and data are in excel- 
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lent agreement at galactocentric distances less than ~25 kpc. However, at larger distances 
from the Galactic center, and outside the Sgr dSph tidal stream, the model overpredicts 
the observed counts by about a factor of two. This mismatch strongly suggests that the 
halo stellar nu mber densit y profi le is much steeper beyond 30 kpc than given by the best- 
fit model from I Juric et al.l (120081 ) (which was constrained using stars at distances below 20 
kpc). Unfortunately, due to the small stripe 82 area, and the effect of the Sgr dSph tidal 
stream on measured stellar counts, it is not possible to robustly derive precise quantitative 
adjustments to model parameters (e.g., there is a strong degeneracy between oblateness and 
density profile power-law index over such a small area). Nevertheless, this analysis provides 
the first direct evidenc e that the ha l o stel lar number density is about a factor of two smaller 
than predicted by the lJuric et al.l (120081 ) model at ~30 kpc from the Galactic center. As 
a rough estimate, in order to induce a factor of two difference in the count ratio as the 
galactocentric radius increases from 25 pc to 35 kpc, the power-law index has to decrease by 
about two (e.g., from to R~^). This estimate is consistent with the profiles shown in 
Figure [H] (beyond 30 kpc). 

The bottom right panel in Figure \2M shows another overdensity at 310° <R.A.< 330° 
and distances in the range 10-25 kpc (about a factor of 1.6 compared to the smooth halo 
model; the Sgr dSph tidal stream produces a peak overdensity of about a factor of five). This 
overdensity appears to be the same structure as the Hercules-Aquila Cloud discovered by 
Belokurov et al.l (j2007bl ). and also seen in the distribution of RR Lyrae stars (Section l4.3p . 



5.1. The Metallicity of the Sgr dSph Tidal Stream (Trailing Arm) 



Ivezic et al.l (120081 ) developed a photometric metallicity estimator based on the u — g 
and g — r colors of main sequence stars. The deep co-added stripe 82 data can be used 
to apply this method to a larger distance limit (~20-30 kpc) than for nominal SDSS data 
(~10 kpc), due to the decreased photometric errors. Figure [251 shows the median metallicity 
as a function of apparent magnitude and position along the stripe 82. Throughout most 
of the volume, the median metallicity is [Fe/H] ~ —1.5, in agreement with results from 



Ivezic et al.l (120081 ). However, in the region where the highest concentration of Sgr dSph 
tidal stream stars is found, the median metallicity is higher by 0.3 dex. On the other hand, 
the median metallicity in the Hercules-Aquila Cloud region is consistent with the median 
halo metallicity. 

In regions outside the Sgr tidal stream, and for r > 21, there is evidence for a metal- 
licity gradient, and for median values [Fe/H] < —1.5. We have analyzed the metallicity 
distribution for stars with 320° <R.A. < 350° (/ ~ 60°, b 40°, 0.3 < g - i < 0.4), 



-27- 



as a function of galactocentric radius. In the range 10-20 kpc, we find that the median 
metaUicity decreases nearly linearly from [Fe/H] = —1.4 (with a rms scatter of 0.33 dex) 
at 10 kpc, to [Fe/H] = —1.6 (with a rms scatter of 0.8 dex, dominated by measurement 
errors) at 20 kpc. T his is a larger gra dient (0.02 dex/kpc) than the upper limit of 0.01 



dex/kpc measured by llvezic et al.l (120081 ) at smaller galactocentric distances. This change of 
metaUicity corresponds to a change of median u — g color (which provides the signal) from 
0.92 mag at 10 kpc (rms scatter around the median is 0.06 mag) to 0.88 mag at 20 kpc (rms 
scatter around the median is 0.17 mag). Clearly, this measurement would greatly benefit 
from deeper w-band data. 

The photometric metaUicity measurement for the Sgr dSph tidal stream can be used 
to test a novel method for estimating the metaUicity of spatially coherent substructures 
when the w-band data are not available. The absolute magnitudes of RR Lyrae and main 
sequence stars have opposite dependence on metaUicity, with RR Lyrae getting brighter and 
main sequence stars getting fainter with decreasing metaUicity (see Equation [TT] in this work 



and Equation A2 in llvezic et al.l (120081 ). respectively). With an assumption that both RR 
Lyrae stars and main sequence stars are at the same distance and have the same metaUicity 
in a well-defined overdensity such as the Sgr dSph tidal stream, then both distance and 
metaUicity can be simultaneously determined. The application of this method to the Sgr 
dSph tidal stream (trailing arm) data is illustrated in Figure | 



A slight difficulty in the application of this method comes from the small RR Lyrae 
sample size and the finite width of their apparent magnitude distribution. We have adopted 
the mode (with 10 stars in a 0.5 mag wide bin) as a relevant statistic. The main sequence 
stars are so numerous that the median and mode are indistinguishable. A consistent solution 
is obtained for a distance modulus of 17.2 mag (distance of 27.5 kpc) and [Fe/H] = —1.2. 

It is encouraging that this metaUicity is the same as obtained using the photometric 
metaUicity method. The photometric metaUicity estimate implies that, with the various 
error distributions taken into account, the metaUicity difference between RR Lyrae and 
main sequence stars in the Sgr trailing arm is not larger than 0.2 dex. We note that the 
photometric metaUicity scale is tied to the SDSS spectroscopic metaUicity scale for main 
sequence stars. 

In addition to similarity of their metaUicity distributions, RR Lyrae and main sequence 
stars in the Sgr trailing arm seem to have similar age distributions. The top middle panel in 
Figure 1251 robustly shows the main sequence turn-off (MSTO) of the trailing arm at q—i ~ 0.3 



(F dw arfs). Together with the color- metallicity-age map shown in Figure 28 from lAn et al. 



(12009! ) . this MSTO color and our metaUicity estimate {[Fe/H] = —1.2) indicate that the 



main sequence population in the Sgr trailing arm is at least 8 Gyr old, and thus consistent 



- 28 - 



with the presumed age of RR Lyrae stars. 

In summary, if both RR Lyrae and main sequence stars are detected in the same spatially 
coherent structure, the metallicity can be reliably estimated (to within ~0.1 dex) by requiring 
the same distance modulus. We note that the calibration of both relevant absolute magnitude 
vs. metallicity relations is based on available globular cluster data. 



Conclusions and Discussion 



We have presented the most complete sample of RR Lyrae stars identified so far in the 
SDSS stripe 82 data set, consisting of 379 RRab and 104 RRc stars. Our visual inspection 
of single-band and color light curves insures that the sample contamination is essentially 
ne gligible. Although the s ky area is re l atively small compared to other recent surveys, such 



as 



Keller et al.l (120081 ) and lMiceli et al.l (120081 ). this RR Lyrae sample has the largest distance 
limit to date (~100 kpc). 

In addition, a large number of well-sampled light curves in the ugriz photometric system 
have enabled the construction of a new set of empirical templates. This multi-band empirical 
template set provides strong constraints for models of stellar pulsations, through the single- 
band light curve shape, the shape variation with wavelength, and through the distributions 
of the templates. For example, while a6-type RR Lyrae appear to have a continuous dis- 
tribution of light curve shapes, c-type RR Lyrae exhibit a clear bimodal distribution. This 
bimodality suggests that the latter class may include two distinct subclasses, which would 
be an unexp ected result g iven that RR Lyrae stars have been carefully studied for over a 
century (e.g., lBaileylll903l ). A subclass of RR Lyrae stars with amplitudes and periods simi- 
lar to the c-type, bu t with asymmetric light curves, has been noted in the MACHO data by 
Alcock et al.l (119961 ). They suggested that these stars are pulsating in the second-harmonic 
mode. 

Our accurate templates will greatly assist the identification of RR Lyrae stars in data sets 
with relatively small numbers of epochs. For example, all three major upcoming wide-area 
ground-based surveys plan to obtain fewer than 10 epochs (SkyMapper, Dark Energy Survey 
and Pan-STARRSl Svr survey), in the same photometric system. Without template fitting, 
the selection based on only low-order light-curve statistics includes a significant fraction 
(~30%) of contaminants (dominated by S Scuti stars). 

The high level of completeness and low contamination of our resulting sample, as well as 
the more precise distance estimates, enabled a more robust study of halo substructure than 
was possible in our first study. Our main result remains: the spatial distribution of halo 
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RR Lyrae at galactocentric distances 5-100 kpc is highly inhomogeneous. The distribution 
of p/Pmodei sug gests that at leas t 20% of the halo within 30 kpc is found in apparently real 
substructures. iBell et al.l (120081 ) find a lower limit of 40%, but they used a different metric 
and different tracers (turn-off stars). Schlaufman et al. (in prep) have demonstrated that at 
least 34% of inner-halo (in their paper they explored regions up to 17 kpc from the Galactic 
center) main sequence turnoff stars are contained in structures they refer to as ECHOS 
(elements of cold halo substructures), based on an examination of radial velocities in the 
SEGUE sample. Thus, three separate studies, using very different probes, have come to 
essentially identical conclusions. 

A comp arison of the observ ed spatial distribution of RR Lyrae stars and main sequence 
stars to the lJuric et al.l (120081 ) model, which was constrained by main sequence stars at 
distances up to 20 kpc, strongly suggests that the halo stellar number density profile steepens 
beyond ~30 kpc. While various indirect evidence for this behavior, based on kiriemati cs of 
field stars, globular clusters, and other tracers has been published ( ICaroUo et al.l 120071 : and 
references therein), our samples provide a direct measurement of the stellar halo spatial 
profile beyond the galactocentric distance limit of ~30 kpc. A s imilar steepe r iing o f the 
spatial profile was detected using candidate RR Lyrae stars by iKeller et al.l (120081 ). In 
addition to presenting new evidence for the steep ening of halo stella r number density profile 
beyond ~30 kpc, we have confirmed the result of iMiceli et al.l (120081 ) that the density profile 
for Oosterhof II stars is steeper than for Oosterhof I stars within 30 kpc. 

We have analyzed several methods for estimating the metallicities of RR Lyrae stars. 
Using spectroscopic data processed with the SDSS SSPP pipeline as a training sample, we 
have established a weak correlation with the u — g color at minimum light. The rms scatter 
around the best-fit relation is ~ 0.3 dex, but systematic errors could be as large, and the 
best-fit relation should be used with care. We have established a relationship between the 
parametrization of our template set for a6-type RR Lyrae stars and the results based on the 
Fourier expansion of light curves. This allowed us to estimate the metallicity from observed 
light curves with an estimated rms scatter of ~ 0.3 dex. This method failed to uncover 
a systematic metallicity difference between field RR Lyrae stars and those associated with 
the Sgr stream, which is seen in main sequence stars. Nevertheless, it is likely that this is 
a promising method for estimating metallicity, which is currently limited by the lack of a 
large and reliable calibration sample. Our templates can be used to bypass noisy F ourier 
tran sformation of sparsely sampled data, as already discussed by iKinemuchi et al.l (120061 ) 
and lKovacs fc Kupil J2007h . 



We introduced a novel method for estimating metallicity that is based on the absolute 
magnitude vs. metallicity relations for RR Lyrae stars and main sequence stars (calibrated 
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using globular clusters). This method does not require the w-band photometry, and will be 
useful to estimate metallicity of spatially coherent structures that may be discovered by the 
Dark Energy and Pan-STARRS surveys. While the existing SDSS data are too shallow to 
apply this method to the Pisces stream, we used it to obtain a metallicity estimate for the Sgr 
tidal stream that is consistent with an independent estimate based on the photometric w-band 
method for main sequence stars. Our result, [Fe/ H] = —1.2, with an uncertainty of ~0.1 dex, 
strongly rules out the hypothesis that the trailing arm of the Sgr dSp h tidal stream has th e 
same metalli c ity as halo field stars {[Fe/H] = —1.5), as suggested by lWatkins et al.l (120091 ). 
Yannv et al.l J2009h detected peaks at [Fe/H] = -1.3 and at [Fe/H] = -1.6 using SDSS 
spectroscopic metallicities for blue horizontal-branch (BHB) stars from the trailing arm. 
However, both the Yanny et al. and Watkins et al. results could be affected by systematic 
errors in the SDSS metallicity scale for BHB stars. Our results suggest that the SDSS 
metallicity scale for BHB stars could be biased low by about 0.3 dex (relative to the SDSS 
metallicity scale for main sequence stars, and assuming that RR Lyrae and main sequence 
stars from the Sgr tidal stream in stripe 82 area have the same metallicity distributions). 



Simulations by lBuUock fc JohnstonI (l2005l ) and I Johnston et al.l (l2008l ) predict that there 
should be a difference in the chemical composition between stars in the inner halo that 
was built from accretion of massive satellites about 9 Gyr ago, and outer halo dominated 
by stars coming from dS ph satellites that were accreted in the last 5 Gyr (see Fig. 11 in 
Bullock fc Johnstorul2005l ). These accreted dSph satellites were pres umably more metal-poo r 
than the massive satellites accreted in earlier epochs (see Fig. 3 in [Robertson et al.l 120051). 
Other recent simulation studies support these conclusions. For example, iDe Lucia &: Helmi 
( 120081 ) find evidence in their simulation for a strong concentration of (relatively) higher 
metalhcity stars at distances close to the Galactic center, and t he presence of ( relati vely) 
lower metallicity stars at distances beyond 20 kpc from the center. IZolotov et al.l (12009! ) find 
from their simulations that their inner halos include stars from both in-situ formed stars and 
accreted populations, while their outer halos appear to originate through pure accretion and 
disruption of satellites. Salvadori et al. (in prep) have considered the distribution of ages 
and metallicities of metal-poor stars in a Milky- Way like halo, as a function of galactocentric 
radius, based on a hybrid N-body and semi-analytic simulation. They find an inner-halo 
population that is well- described by a power-law index n = —3.2 (for stars with —2 < 
[Fe/H] < —1), and an outer-halo consistent with a much shallower profile, n = —2.2. The 
relative contributions of stars with [Fe/H] < —2 in their simulation increases from about 
16% for stars within 7 kpc < i? < 20 kpc, to > 40% for stars with R > 20 kpc. 

Our observational results provide some support for these simulation-based predictions. 
The faintest main sequence stars (r ~ 21.5), not including those apparently associated 
with the Sgr stream, exhibit median metallicities at least 0.3 dex lower {[Fe/H] = —1.8, 
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see the left panel in Figure 1251) than halo stars within 10-15 kpc from the G alactic center 



(Ivezic et al. 



20081 ). Other studies provide additional support. For example, ICaroUo et al. 



( 120071 ) and CaroUo et al. (in prep) have indicated that field stars likely to be associated with 
the outer halo exhibit metallicities that are substantially lower {[Fe/H] = —2.2) than those 
of the inner halo (they place the inner/outer halo boundary at ~10 kpc). The two metallicity 
measurements are fully consistent because photometric w-band metallicity estimator is biased 
high for spectroscopic values with [Fe/H] < —2 (Bond et al., in prep). 

Carollo et al. (in prep) also found that the inferred density profiles of the inner- and 
outer-halo populations differ as well; the inner halo being consistent with a power-law profile 
with index n ~ —3.4, and the outer halo having index n = —2.1 (the data analysed by Carollo 
et al. are not suitable for determination of the relative normalizations of the inner- and 
outer-halo populations, due to the selection effects involved). This difference in the density 
profiles is very similar to the difference in density profiles for Oosterhof I and Oosterhof II 
stars. However, we emphasize that the steeper profile for the inner halo from Carollo et 
al. corresponds to more metal-rich stars, while we found weak evidence that Oosterhof II 
stars tend to be more metal poor than Oosterhof I stars. 

Simulations indicate that high surface brightness substructures in the halo originate 
from single satellites, typically massive dSph which tend to be accreted over the last few 
Gyr (IBuUock fc Johnston! |2005| ) . and these massiy e galaxies are expect ed to be mor e more 



metal-rich than halo field stars (IFont et al.l l2008l ) . The results from llvezic et al.l (120081 ) 



and the results presented here seem to support this prediction. The inner halo has a median 
metallicity of [Fe/H] = —1.5, while at least two strong overdensities have higher metallicities 
- the Monoceros stream has [Fe/H] = —1.0, and for the trailing part of the Sgr tidal stream 
we find [Fe/H] = —1.2. We emphasize that these three measurements are obtained using 
the same method/calibration and the same data set, an d thus the measureme nts of relative 
differences are expected to be robust. Simulations by [Johnston et al.l (120081 ) predict that 



the m etallicity of low surface brightness features, such as the Virgo Overdensity (jJuric et al. 



20081). is expected to be lower than the median halo metallicity. Although the initial estimate 
by lJuric et al. J2008h ( see their Figure 39) claimed a metallicity for this structure of [ Fe/H] 



— 1.5, a recent study of the photometric metallicity of stars in this structure by lAn et al. 
( 120091 ) have suggested that the mean metallicity is even lower: [Fe/H] = —2.0. Both of 



these studies are limited by potentially large systematics associated with the photometric 
metallicity technique, so more detailed analysis of spectroscopically determined metallicities 
for stars in this structure would be of great value. The sparse information that is available, 
based on spectroscopic determinations, supports a mean metallicity between [Fe/H] = — 1. 8 
and -2.2 JPuffau et a l."20od: Vivas et ahlbood : IPrior et al.lbood : ICasetti-Dinescu et al.l[2009[ ) . 
Casetti-Dinescu et al.l (i2009l ) report the measurement of a precise absolute proper motion 
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for the RR Lyrae star RR 167, which appears highly hkely to be associated with the Virgo 
Overdensity. This proper motion, in combination with their distance estimate (17 kpc from 
the Galactic center) and radial velocity measurement, indicate that the Virgo Overdensity 
structure may well be on a highly destructive orbit, with pericenter ~ 11 kpc, and apocenter 
~ 90 kpc. Thus, the interpretation of the Virgo Overdensity as a dwarf galaxy, perhaps 
similar to those that participated in formation of the outer-halo population, is strengthened. 

Our result that (inner-) halo stellar number density profile steepens beyond ~30 kpc is 
limited by the relatively small distance limit for main sequence stars (35 kpc), the sparseness 
of the RR Lyrae sample (~500 objects), and the small survey area (~300 deg^). Ideally, the 
halo stellar number density profile should be studied using numerous main sequence stars 
detected over a large fraction of sky. To do so to a distance limit of 100 kpc, imaging in 
at least g and r bands (or their equivalent) to a depth several magnitudes fainter than the 
co-added SDSS stripe 82 data is required (r > 25). Pan-STARRS, the Dark Energy Survey 
and LSST are planning to obtain such data over large areas of sky. The LSST, with its deep 
M-band data, will also e xtend metalli c ity m apping of field main sequence stars over half of 



the sky in the south; see llvezic et al.l (l2008l ) for details. For substructures to be potentially 
discovered in the north by Pan-STARRS, the method presented here can be used to estimate 
the metallicity of spatially coherent structures even without the w-band data. 
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Table 1. Positions and ugriz photometry for 483 stripe 82 RR Lyrae stars 





Dec°- 


"■MJD 




^err 


9m JD 


9" 


^err 


' MJD 




err 


'•MJD 




err 


^MJD 


2^= 


'^err 


0.935679 


1.115859 


51075.3024 


18.04 


0.01 


51075.3041 


16.79 


0.01 


51075.3007 


16.65 


0.01 


51075.3016 


16.62 


0.01 


51075.3032 


16.59 


0.01 


0.935679 


1.115859 


52196.3030 


18.54 


0.02 


52196.3047 


17.35 


0.01 


52196.3013 


17.05 


0.01 


52196.3022 


16.94 


0.01 


52196.3038 


16.88 


0.01 


0.935679 


1.115859 


52197.3161 


-99.99 


-99.00 


52197.3178 


17.21 


0.01 


52197.3144 


16.89 


0.01 


52197.3153 


16.77 


0.01 


52197.3169 


16.75 


0.01 


0.935679 


1.115859 


52225.1663 


18.53 


0.02 


52225.1679 


17.34 


0.01 


52225.1646 


17.05 


0.01 


52225.1654 


16.94 


0.01 


52225.1671 


16.90 


0.01 


0.935679 


1.115859 


52231.2020 


-99.99 


-99.00 


52231.2037 


16.98 


0.01 


52231.2003 


16.77 


0.01 


52231.2012 


16.67 


0.01 


52231.2029 


16.65 


0.01 



^Median equatorial J2000.0 right ascension and declination in degrees 
''Modified Julian date in days 

■^Magnitude is not corrected for ISM extinction, -99.00 if the measurement is unreliable 
•^Magnitude uncertainty, -99.00 if the measurement is unreliable 



Note. — Table[T]is published in its entirety in the electronic edition of the Journal. A portion is shown here for guidance regarding its form and content. 



Table 2. Light curve parameters for stripe 82 RR Lyrae stars 



RR Lyrac^ 


Type 


Period'' 


Al 


"0 


't'o'" 


Ti 








^/ 




^0 




t/ 


Af 


'o 


't'o" 


t/ 


Al 


^0 


't'o' 


Tf 


4099 


ab 


0.642 


0.53 


17.82 


51075.27 


101 


0.56 


16.66 


51075.28 


101 


0.40 


16.56 


51075.29 


103 


0.32 


16.55 


51075.29 


102 


0.30 


16.54 


51075.28 


100 


13350 


ab 


0.548 


1.03 


18.24 


54025.33 


107 


1.09 


16.99 


54025.32 


120 


0.80 


17.09 


54025.32 


112 


0.64 


17.15 


54025.32 


116 


0.58 


17.19 


54025.32 


114 


15927 


ab 


0.612 


0.65 


18.92 


53680.22 


102 


0.70 


17.65 


53680.22 


106 


0.49 


17.60 


53680.22 


108 


0.37 


17.61 


53680.24 


104 


0.35 


17.62 


53680.24 


100 


20406 


ab 


0.632 


0.55 


16.38 


54000.28 


100 


0.59 


15.21 


54000.29 


101 


0.43 


15.12 


54000.27 


108 


0.34 


15.12 


54000.29 


102 


0.30 


15.13 


54000.29 


100 


21992 


ab 


0.626 


1.11 


15.53 


53698.24 


106 


1.14 


14.38 


53698.24 


112 


0.85 


14.42 


53698.24 


114 


0.66 


14.52 


53698.24 


111 


0.62 


14.52 


53698.24 


114 



^RR Lyrac ID number 
'^Period in days 

'^Amplitude determined from the best-fit template 

'^Mean magnitude corrected for the ISM extinction and determined from the best-fit template 
'^Epoch of maximum brightness determined from the best-fit template 
^Best-fit template ID number 

Note. — Table [2] is published in its entirety and with higher-precision values in the electronic edition of the Journal. A portion is shown here for guidance regarding its form and content. 



00 

00 



Table 3. Positions and derived parameters of stripe 82 RR Lyrae stars 



RR Lyrae" 


R.A.'' 


Dec'' 


rExt" 




"-GC 












{vr 


(u ~ q)^ - 




(q — r)^ - 


{S-r)min 


4099 


0.935679 


1.115859 


0.089 


17.75 


20.03 


18.13 


16.99 


16.78 


16.70 


16.68 


16.85 


1.10 


0.02 


0.28 


0.01 


13350 


0.283437 


1.178522 


0.080 


24.77 


26.55 


18.84 


17.68 


17.54 


17.50 


17.50 


17.57 


1.16 


0.01 


0.26 


0.01 


15927 


3.254658 


-0.584066 


0.090 


29.12 


30.96 


19.29 


18.06 


17.86 


17.79 


17.78 


17.92 


1.24 


0.02 


0.28 


0.01 


20406 


3.244369 


0.218891 


0.088 


9.13 


12.76 


16.71 


15.54 


15.34 


15.29 


15.28 


15.40 


1.18 


0.01 


0.27 


0.01 


21992 


4.315354 


1.054582 


0.077 


7.35 


11.54 


16.19 


15.04 


14.91 


14.86 


14.85 


14.93 


1.13 


0.01 


0.25 


0.01 



''RR Lyrae ID number 

'^Median equatorial J2000.0 right ascension and declination in degrees 

'^ISM extinction in the SDSS r band obtained from lSchlegel. Finkbeiner fc David lll998l ') 

'^Heliocentric and galactocentric distance in kpc 

°Flux-averaged magnitudes corrected for ISM extinction 

^Median u — g color (corrected for ISM extinction) at the minimum brightness and the uncertainty in the median 

s Median g — r color (corrected for ISM extinction) at the minimum brightness and the uncertainty in the median ^ 



Note. — 



Tabic [3] is published in its entirety in the electronic edition of the Journal. A portion is shown here for guidance regarding its form and content. 
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Table 4. ISesar et al.l (120071 ) overdensities revisited 



S07 label" 




NRRLyr 


M<1 


Are 

non — Liar 


NRRLyr/Ni^^ 


A 


84 


76 


7 


1 


0.91 


B 


144 


98 


44 


2 


0.68 


C 


54 


50 


4 





0.93 


D 


8 


5 


3 





0.63 


E 


11 


7 


4 





0.64 


F 


11 


4 


7 





0.36 


G 


10 


5 


5 





0.50 


H 


7 





7 





0.00 


I 


4 


1 


3 





0.25 


J 


26 


12 


14 





0.46 


K 


8 


1 


7 





0.13 


L 


3 





3 





0.00 


M 


5 


1 


4 





0.20 


all 


375 


260 


112 


3 


0.69 



^Overdensity label from S07 

'^Number of S07 candidate RR Lyrae in the overdensity 
'^Number of true RR, Lyrae stars 
'^Number of variable, non-RR Lyrae stars 
''Number of non- variable stars 
'Fraction of true RR Lyrae stars 
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Table 5. Poisson statistics for RR Lyrae detection in six stripe 82 regions 



Region 




oos 


Statistical significance'* 


30° < R.A. < 60° 


12 ±4 


2 


2.5 


50 < d/kpc < 75 0° < R.A. < 30° 


10 ±3 


6 


1 


310° < R.A. < 340° 


19 ±4 


8 


3 


30° < R.A. < 60° 


9±3 





3 


75 < dkpc < 100 0° < R.A. < 30° 


8± 3 


3 


1.5 


310° < R.A. < 340° 


14 ±4 


4 


2.5 



^Expected number of RR Lyrae calculated by integrating Eq. 1161 
''Uncertainty in the number of expected RR Lyrae calculated as ^ Ne: 
Observed number of RR Lyrae 

Statistical significance in units of a, calculated as \Nc,xp — N^haM'^ 
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Table 6. Positions and derived parameters of RR Lyrae stars in the Sagittarius and Pisces 

streams 



Stream 


{R.A.)<^ 




ao 


N- 


Ool 


A7"c 

"OoII 






Sagittarius 


30° 


25 


38 


7 


29 


9 


0.6 


1.13 ±0.01 


Pisces 


355° 


80 


15 


1 


15 





0.6 


1.12 ±0.06 



''Median equatorial J2000.0 right ascension 

''Median heliocentric distance in kpc 

'^Number of type ab, c, Oo I, and Oo II RR Lyrae stars 

"^Median period of RRab stars in days 

''Median {u — g)min color of RR Lyrae stars and the uncertainty in the median 



1 10 100 

time between succesive visits (days) 

Fig. 1. — The cadence of SDSS stripe 82 observations. On average, the objects were most 
often re-observed every two days (the SDSS-II SN Survey cadence), followed by 5-day, 10-day 
and yearly re-observations. 
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Fig. 2. — The distribution of (7-band amplitudes (Ag) and mean magnitudes (qq) determined 
from best-fit templates for 3449 RR Lyrae candidates. The solid line separates candidates 
into template-accepted {blue dots) and ambiguous candidates {green triangles), which lie 
above the line, and template-rejected {red open squares) which lie below the line. 
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Fig. 3. — Top: The ag vs. distribution of template-accepted {blue dots), template- rejected 
{red open squares), and ambiguous candidates {green triangles) before the \ag — lA2ar\ < 0.03 
selection (region between dashed lines). Note that the majority of template-rejected and 
ambiguous candidates do not follow the cfg = 1.420"^ relation typical of RR Lyrae stars {solid 
line), but rather the ag = ar relation {dotted line) that indicates non-pulsational variability. 
Middle: The Cg-r vs. sqrt{(jg + af.) distribution of candidates from the top panel after the 
\(Tg — 1.42(Tr.| < 0.03 selection is applied. The dashed lines show the selection based on the 
ag-r and sqrt{ag + a^) values, where ag^r is the rms scatter in the g — r color. Bottom: 
The Ag vs. {g) distribution of candidates selected from the middle panel's dashed box region. 
Note that almost all of the ambiguous candidates with g^ ~ 20.5 shown in Fig.[2]are removed 
after top and middle panel cuts are made, while the template-accepted candidates remain. 
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19.6 E 




phase ph 



Fig. 4. — The top and middle panel show the g- {left) and r-band {right) light curves of two 
RR Lyrae stars with the observed ag slightly smaller than the predicted ag = 1.42(7^ value 
(by ~ 0.1 mag). The bottom panel shows the g- and r-band light curves of a RR Lyrae star 
with ag-r greater by ~ 0.05 mag than the cut-off dashed hne in Figure [3] ( bottom) . The solid 
lines show the best-fit g- and r-band templates (see Section [3]). All three stars have slightly 
peculiar ag, ar, and cTg^^r values, but are still included in the final RR Lyrae sample as their 
light curves are very much RR Lyrae-like. 
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RR Lyrae 3094190, P = 0.21, A = 0.37, template = 8 




phase 



phase 



Fig. 5. — The g {top) and g — r [bottom) light curves of a RRc star [left) and a candidate 
echpsing binary (right). The best fit V-band templates are shown as solid lines. Even though 
their g-hand light curves are quite similar, and their periods are in the 0.2 < P < 0.43 days 
range typical of RRc stars, their g — r light curves are quite different. The eclipsing binary 
is outside the Figure [3] (middle) selection box as it has a much smaller rms scatter in the 
g — r color (dg-r ~ 0.02 mag) than the RRc star ((Jg^r ~ 0.07 mag). 
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phase08 : 0.824 to 0.911 




Fig. 6. — Top: The 22 ^f-band RR a6-type template light curves, color-coded from black to 
red using 0q g, defined as the phase at which a q-b and template value is equal t o O.8. The 
(t)n.H is s i milar in principle to the rise time used by ISandage. Katem fc Sandagd (1l98ll ) and 
LaydenI (119981 ) (shorter rise time means higher ^o.s)- Bottom: A comparison of two g band 



RR c-type template light curves. The template showed by the dashed line has a higher 
value than the template showed by the solid line (0.80 vs. 0.72). 



0.8 
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0.0 



0.2 



0.4 0.6 
phase 



0.8 



.0 



Fig. 7. — An example of incompleteness in our RR Lyrae light curve template set. The 
period-folded g- and r-band light curves are shown with symbols, and the best-fit templates 
are shown with solid lines. The fact that even the best-fit templates fail to model these light 
curves at minimum brightness (0.5 < phase < 0.9) suggests that our light curve template 
set is not entirely complete. 



0.0 0.2 



0.4 0.6 

phase 



0.8 1.0 



0.0 0.2 



0.4 0.6 

phase 



0.8 1.0 




0.4 0.6 

phase 



Fig. 8. — A comparison of empirical {blue lines) and theoretical ( red lines) RR a 6-tvpe 
template light curves in SDSS ugriz band-passes. The theoretical iMarconi et al.l (120061 ) 
templates were obtained from a RR Lyrae pulsational model with metallicity Z = 0.001, 
luminosity L = I.6IL0, mass M = 0.75Mq, and effective temperature Tg// = 6100 — 6900 
K. Note how the light curve shape of empirical and theoretical templates changes significantly 
between griz bands, becoming less convex (less curved upward) on the descending brightness 
branch (phase< 0.7) as the band-pass moves to longer wavelengths. The shape does not 
change significantly between u and g bands. 
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CP 



0.3 

0.2 

0.1 

0.0 

0.1 
0.2 
0.3 

0.2 

0.1 

0.0 

0.1 
0.2 
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phase 







0.5 



. 1 I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I LJ I I I I I I I I I I I I 




I I ipi I I I I I I I I I I I I I I I I 




0.9 1.0 1.1 1.2 0.9 1.0 1.1 1.2 

u — g u — g 



Fig. 9. — The u — g vs. g — r color-color loops of four RR Lyrae stars. The top left panel 
shows a long-period, low-amplitude type b RR Lyrae stars (ID 276162), the top right panel 
shows a short-period, large-amphtude type a RR Lyrae star (ID 293282), and the bottom 
panels show two type c RR Lyrae stars (IDs 376465 and 429508), with the one on the left 
having a longer period and a lower amplitude than the one on the right. The observed 
{symbols with error bars) and predicted u — g and g — r colors (from the best-fit ugr-ha.nd 
templates, lines) are color-coded according to the phase of pulsation. The error bars are 
not shown in three panels to avoid cluttering. The amplitude in the g — r color increases 



with the (7-band amplitude {Ag 



OAAg, since Ag ~ 



1.4 A, and An 



Ag — Ar 



The 



color-color loops for the rest of the RR Lyrae stars are provided in the electronic edition of 
the Journal. 
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0° 




o 



Fig. 10. — The top panel shows the spatial distribution of 366 RRab stars from stripe 82, 
while the bottom panel shows 302 stars randomly drawn from a smooth model distribution 
(Eq. [T6|) . Both samples of stars have heliocentric distances between 5 and 120 kpc and 
\Dec\ < 1.25°. The radial axis is the heliocentric distance in kpc and the angle is the 
equatorial right ascension. The circles correspond to {V) of 16.3, 17.8, 18.7, and 19.3 mag 
(corrected for ISM extinction). The RRab stars shown in the top panel are divided into 
Oosterhof I (289 stars, red triangles) and Oosterhof II (77 stars, blue dots) using the selection 
boundary shown in Fig. [161 








Fig. 11. — The number density distribution of observed (top) and smooth model-based 
{bottom) RR Lyrae stars from Fig. [TOl computed using Eq. [181 shown relative to the 
smooth model number density (Pmodei)- The log(p/ p^dei) values are color-coded according to 
the legend, with values outside the legend's range saturating in purple and red, respectively. 
The stripe 82 plane intersects the Sagittarius dSph tidal stream (trailing arm) at (~ 25 
kpc, ~ 30°), the Hercules-Aquila cloud at (< 25 kpc, ~ 330°), and the Pisces stream at 
(~ 80 kpc, ~ 355°). The overdensities detected in S07 are overplotted and labeled in the 
top panel. The "+" symbol marks the position where the lLaw. Johnston fc Majewskil (120051 ) 
"spherical" model of the Sagittarius dSph tidal stream (trailing arm) intersects the stripe 82 
plane. 



Fig. 12. — A single frame of an animation showing a fly- by of the stripe 82 plane 
and of the scaled Galactic plane (annotated artist's concept by NASA/JPL-Caltech). 
The white dots show the Sagitt arius dSph and its tidal streams, as modeled by the 
Law. Johnston fc Majewskil (120051 ) "spherical" model. The animation is provided in the 



electronic edition of the Journal. 
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Fig. 13. — The differential {top, log y-axis) and cumulative {bottom) distribution of 
log(p/p^^^g;) in model-based {dashed line) and observed number density maps {solid line, 
within 110 kpc; dotted line, within 30 kpc). The slope of the observed log(p/p^^^g;) distribu- 
tion in the 0.3 < log(p/p^^^g,) < 1.0 range is -1.65 ±0.05 within 110 kpc, and -1.91 ±0.09 
within 30 kpc. Note significantly higher fractions of strongly overdense and underdense 
regions than predicted by a smooth halo model. 
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Fig. 14. — The radial number density profiles for four regions in R.A. : 30° < R.A. < 60° 
{top left), 0° < R.A. < 30° {top right), 340° < R.A. < 0° {bottom left), and 310° < R.A. < 
340° {bottom right). The bottom left panel contains the Pisces overdensity at Rgc ~ 80 
kpc (log(i?Gc) ~ 1-9). The symbols show the observed median number density and the 
rms scatter around the median value for pixels shown in Figure [11] {top). The dashed 
line shows the number density predicted from the smooth model (log(p^^^g;),Eq. [16]), and 
the dotted lines show the ±0.3 dex envelope (a factor of 2) around the predicted values 
which corresponds to the sample Poisson noise de termined from naodel- based number density 
maps. The sohd line shows the prediction from IWatkins et al.l (12009). The p^,^ in Eg. [16] 
was increased by 33% to 5.6 kpc~^ to match the observations, and [Watkins et al.l (120091 ) 
predictions (their Eq. 12) had to be scaled down by a factor of 10 (assuming that their 
normalization was meant to be in units of stars per kpc^). The vertical line is added to 
guide the eye and shows the Rqc = 30 kpc. 
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Fig. 15. — Top: The radial number density profiles of Oo I {crosses) and Oo II RRab stars 
{dots) in the 340° < R.A. < 30° direction. The symbols show the observed median number 
density and the rms scatter around the median value. The best-fit po wer laws to the data 
(power law indices of —2.3 and —4.5) are shown as solid lines, with the iMiceli et al.l (120081 ) 
power laws overplotted as dashed lines. Bottom: The radial number density profiles of RRab 
{dots) and RRc stars {crosses) in the 340° < R.A. < 30° direction. The best-fit power laws 
{solid a.iad dashed line) have consistent power law indices (—2.9). 
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Fig. 16. — The period-amplitude diagram for RR Lyrae stars with high {( < 2.3, circles) 
and low quality ^f-band light curves (C > 2.3, stars), color-coded by the median g — r color 
at minimum brightness {{g — r)min) in the top panel, and by 0Qg in the middle panel. The 
vertical line at log{Period) = —0.36 (10.5 hours) divides the RRc (first overtone, to the 
left) and RRab stars (fundamental mode, to the right). The solid line shows the position 
of the Oo I locus, while the dashed line (offset by 0.03 in the log{Period) direction from 
the Oo I locus line) separates the Oo I (to the left) and Oo II RRab stars (to the right). 
The position of the dashed line was determined from the A\og{Period) histogram [bottom), 
where A \og{Period) is the log{Period) distance (at constant amplitude) from the Oo I locus 
line. The g- and r-band light curves for the star with Ag ~ 1.3 and log{Period) ~ —0.09 
are shown in Fig. [71 
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Fig. 17. — The symbols show high quahty (7-band hght curves (normahzed and period- 
folded) of RRc stars fitted with template (eight stars, dots) and 1 (nine stars, triangles) , 
with templates shown as solid lines. The template and 1 are also shown as the solid and 
dashed line in FiguieM (bottom) . This bimodal distribution of observed light curves suggests 
that c-type RR Lyrae stars may include two distinct subtypes. 
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Fig. 18. — The period- amplitude diagram for RR Lyrae stars selected from the Sagittarius 
(dots) and Pisces streams (stars). The vertical line at log{Period) = —0.36 divides the RRc 
(first overtone, to the left) and RRab stars (fundamental mode, to the right). The solid line 
shows the position of the Oo I locus, while the dashed line (offset by 0.03 in the \og{Period) 
direction from the Oo I locus line) separates the Oo I (to the left) and Oo II RRab stars (to 
the right). 
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Fig. 19. — Top: The {u — g)min vs. {g — r)min diagram for 51 RRab stars with high {dots) 
and low quahty (/-band hght curves (stars), and with spectrosco pic metalhcity measurements. 
The spectroscopic metalhcities were obtained by lDe Led (120081) from SPSS s pectra processed 



Lee et al.ll2008l ). The {u-g)r 



with the SEGUE Spectroscopic Parameters Pipehne (SSPP; 
and {g — r)min {symbols) show median colors at minimum brightness (0.5 < phase < 0.7), 
and the error bars show errors in the median. The symbols are color-coded according to 
spectroscopically determined metallicity: [Fe/H] < —2.0 {red), —1.5 < [Fe/H] < —2.0 
{green), and [Fe/H] > —1.5 {black). Bottom: The distribution of Oo I {black symbols) 
and Oo II {red symbols) RRab stars with high {dots) and low quality g-hand light curves 
{stars) in the [Fe/H] vs. {u — g)min diagram. The Oo II stars are on average more metal- 
poor {[Fe/H] ~ —2.0) and bluer in {u — g)min {{u — g)min ~ 1-1) than the Oo I stars 
{[Fe/H] ~ -1.7). 
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Fig. 20. — The dependence of spectroscopic metallicity on median u — g color at minimum 
brightness {{u — g)min) for RRab {red symbols) and RRc {blue symbols) stars with high {dots) 
and fow quahty g-hand hght curves {stars). The best-fit to median [Fe/H] values in 0.02 
wide {u — g)min bins is [Fe/ H]photo = —46.47 + 72 .36(m — q)mir, — 28.98(n — 5')mm {dashed 
line). The spectroscopic metallicity errors are from lDe Led ((20081), but were not used in the 
fit as they are most probably underestimated (they should be ~ 0.3 dex, De Lee, private 
communication). The rms scatter of 76 individual data points around this line is 0.3 dex, 
and is consistent with the expected intrinsic scatter in spectroscopic metallicity. 
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Fig. 21. — Examples of best-fit Fourier expansion with six sine terms to light curve templates 
for afe-type RR Lyrae stars. The solid line shows the template, and the three dashed lines 
show the first three harmonics (offset for clarity). The sum of all six harmonics is shown by 
the dotted line (barely visible and close to the solid line). These four examples span the full 
range of the minimum phase, ranging from 0.75 in the top left panel to 0.87 in the bottom 
right panel. For fixed values of period and amplitude, the phase of the minimum increases 
as the metallicity decreases. 
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Fig. 22. — Stellar counts in deep coadded data from stripe 82. Each panel shows stars from 
a 0.1 mag wide g — i color bin, as marked on top, with the g — i color increasing from top left 
(0.3-0.4) to bottom right (0.8-0.9). Assuming [Fe/H] = —1.5, the median r-band absolute 
magnitude varies from Mr = 4.4 to Mr = 6.6. The counts are shown on a logarithmic scale, 
increasing from blue to yellow and red (varying normalization for each panel to emphasize 
features). The red vertical features visible at the right edge of each panel are due to disk 
stars at r < 20 and halo stars at fainter magnitudes (including the Hercules- Aquila Cloud). 
The stars in the well-defined overdensity at R.A. ~ 20° — 50° belong to the Sgr dSph tidal 
stream (trailing arm). The variation of apparent magnitude at which the overdensity is 
detected with the g — i color shows that these stars are on the main sequence (see fig. [23]) . 
The distance to the overdensity is in the 25-30 kpc range. 
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Fig. 23. — The top left panel shows the stripe 82 footprint in galactic coordinates. The 
corresponding R.A. along the stripe is marked (Dec~0). The other five panels show the r 
vs. g — i color- magnitude (Hess) diagrams for five 25 deg^ large regions from stripe 82 selected 
by R.A. (the range is listed on top of each panel). The top middle panel corresponds to a 
region intersecting the Sgr dSph tidal stream; note the main sequence turn-off (MSTO) at 
g — i 0.3 (F type dwarfs) and a well-defined sub-giant branch. The MSTO at — i ~ 0.3 
indicates that the main sequence stars in the Sgr trailing arm are at least 8 Gyr old. The 
counts are shown on a logarithmic scale, increasing from blue to yellow and red (with varying 
normalization to emphasize features). The dashed lines are added to guide the eye and to 
show the position of main sequence stars with [Fe/H] = —1.2 at a distance of 30 kpc. Note 
the lack of stars with — z ~ 0.4 and r > 22 in most panels. For the median halo metallicity 
of [Fe/H] = —1.5, r = 22 corresponds to 28.7 kpc. The feature at g — i ^ 0.5 and r ~ 15 — 18 
is due to thick disk stars. 
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Fig. 24. — A comparison of the observed distribution of main sequence s tars in stripe 82 and 
predictions based on a smooth oblate halo model from lJuric et al.l (120081 ). The top left panel 
shows the number density of stars with 0.3 < g — i < 0.4 on a logarithmic scale (increasing 
from blue to yellow and red) in a galactocentric cylindrical coordinate system (the Sun is 
at i? = 8 kpc, Z = 0; all data are at negative Z). The top right panel shows the model 
prediction. The bottom left panel shows the In (data/mo del) map, color-coded according to 
the legend shown in the panel. The overdensity centered on [R ~20 kpc, Z ~20 kpc) is the 
Sgr dSph tidal stream (trailing arm). Note that at large distances from the Galactic center, 
and outside the Sgr dSph tidal stream, the model overpredicts the observed counts by about 
a factor of 2 to 3. At galactocentric distances < 20 kpc, the model and data are in excellent 
agreement. The bottom right panel is analogous to the bottom left panel, except that the 
map is shown in stripe 82 coordinates, with the equatorial position angle (R.A.) marked on 
the top. The Sgr dSph tidal stream is seen at 25° <R.A.< 55°, with the maximum value 
of the overdensity at about five. The Hercules-Aquila cloud is seen at 310° <R.A.< 330°, 
and distances in the range 10-25 kpc, and represents a factor of ~1.6 overdensity. The 
region with 340° <R.A.< 5° and distanc e s in t he range 25-35 kpc represents a factor of ~ 2 
underdensity, relative to the lJuric et al.l (120081 ) smooth halo model. 
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Fig. 25. — The left panel shows the median photometric metallicity for 88,000 stars with 
0.3 < g — i < 0.4 and 18 < r < 22.2 in deep coadded data from stripe 82, as a function 
of apparent magnitude and position along the stripe. Note the relatively high metallicity 

{[Fe/H] 1.20) of the Sgr dSph tidal stream region (R.A.~ 30° - 40° and 20 < r < 21). 

On the other hand, the median metallicity in the region with R.A.< —10° and r > 20, where 
the Hercules-Aquila overdensity is found, is consistent with the median halo metallicity 
{[Fe/H] = —1.50). The right panel compares the metallicity distributions for stars with 
20.5 < r < 21 selected from two regions: 2,200 stars with 30° < R.A. < 40° (Sgr dSph tidal 
stream, right histogram) and 4,500 stars with —40° < R.A. < —10° (control halo sample, 
left histogram). The dashed lines are best-fit Gaussians; both have widths a = 0.4 dex and 
are centered on [Fe/H] = —1.50 (halo, left) and [Fe/H] = —1.20 (Sgr dSph tidal stream, 
right) . 
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Fig. 26. — A summary of the constraints on the distance and metaUicity of Sgr dSph tidal 
stream (traihng arm). The dot-dashed hnes show a constraint obtained from the mode 
of the apparent magnitude distribution of RR Lyrae stars, with ±0.1 mag adopted as the 
uncertainty. Diagonal short-dashed lines are analogous constraints obtained from the median 
apparent magnitude of main sequence stars with 0.4 < g — i < 0.5. The vertical long- 
dashed lines mark the median photometric metallicity for main sequence stars, with ±0.1 dex 
adopted as the uncertainty. All three constraints agree if the distance modulus is DM = 17.2 
mag and [Fe/H] = —1.2. A hypothesis that the Sgr dSph tidal stream (trailing arm) has 
the same metallicity as halo field stars {[Fe/H] = —1.5) is strongly ruled out. 



